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Abstract: The 2d gauge theory on the lattice is equivalent to the twisted Eguchi- 
Kawai model, which we simulated at ranging from 25 to 515. We observe a clear 
large N scaling for the 1- and 2-point function of Wilson loops, as well as the 
2-point function of Polyakov lines. The 2-point functions agree with a universal 
wave function renormalization. The large double scaling limit corresponds to the 
continuum limit of non-commutative gauge theory, so the observed large scaling 
demonstrates the non-perturbative renormalizability of this non-commutative field 
theory. The area law for the Wilson loops holds at small physical area as in com- 
mutative 2d planar gauge theory, but at large areas we find an oscillating behavior 
instead. In that regime the phase of the Wilson loop grows linearly with the area. 
This agrees with the Aharonov-Bohm effect in the presence of a constant magnetic 
field, identified with the inverse non-commutativity parameter. 

Next we investigate the 3d X(j)^ model with two non-commutative coordinates and 
explore its phase diagram. Our results agree with a conjecture by Gubser and Sondhi 
in d = 4, who predicted that the ordered regime splits into a uniform phase and a 
phase dominated by stripe patterns. We further present results for the correlators 
and the dispersion relation. In non-commutative field theory the Lorentz invariance 
is explicitly broken, which leads to a deformation of the dispersion relation. In one 
loop perturbation theory this deformation involves an additional infrared divergent 
term. Our data agree with this perturbative result. 

We also confirm the recent observation by Ambj0rn and Catterall that stripes 
occur even in d = 2, although they imply the spontaneous breaking of the translation 
symmetry. 
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1 Introduction 



The ideas of non-commutative space-time and field theories defined on such a struc- 
ture started akcady in the year 1947. At this time the concept of renormahzation 
was not yet well estabhshed and therefore uhraviolet divergences in quantum field 
theory still caused serious problems. To solve these problems or at least weaken 
them Snyder introduced the quantized space-time [1] (see also [2]). 

The plan was to define quantum field theories on a space-time which is smeared 
out at very small length scales. This means that in addition to Heisenberg's un- 
certainty relation between coordinates and momenta there is a uncertainty between 
different coordinates. 

As in the quantization of the classical phase space, space-time can be quan- 
tized by replacing the usual coordinates by Hermitian operators x^, obeying the 
commutator relation 

= iO^^. (1.1) 

The non-commutativity tensor O^j, is a real-valued antisymmetric dxd matrix and 
d is the space-time dimension. Since the coordinate operators do not commute they 
cannot be diagonalized simultaneously and therefore induce the uncertainty relation 

Ax^Ax,>^\e^,\. (1.2) 

This uncertainty implies a quantum structure of space-time and due to the lack of 
points in space-time it then represents an effective ultraviolet cut-off. 

Much later, in 1996, it was shown by Filk [3] that in field theory on a non- 
commutative plane the divergences of commutative field theory still occur. In ad- 
dition to those divergences the authors of Refs. [4, 5] found in 1999 that there is a 
mixing of ultraviolet and infrared divergences. 

The concept of quantized space-time has not been followed further in the early 
days of quantum field theory, since the renormahzation technique became more and 
more successful. It came up again first in the 80 's, when Connes formulated the 
mathematically rigorous [6] framework of non-commutative geometry. In physical 
theories a non-commutative space-time first appeared in string theory, namely in 
the quantization of open strings [7]. In an constant magnetic background field the 
boundary conditions change and the zero momentum modes of the string do not 
commute anymore. Instead they obey a commutation relation of the type (1.1), 
where |0^i/| is proportional to the inverse background field. 
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The zero momentum modes of an open string can be interpreted as the end points 
of the string, which are confined to a submanifold, i.e. a D-brane. The commutator 
(1.1) imphes a non-commutative geometry on the branes. Hence a quantized space- 
time appears naturally in string theory [8]. 

Another field of interest, where the non-commutativity of space-time could play 
an important role, is quantum gravity. It is an old believe that for a quantized theory 
of gravitation space-time has to change its nature on very small length scales. The 
synthesis of the principles of quantum mechanics and of classical general relativity 
leads to a space-time uncertainty [9,10], which implies that any theory of quantum 
gravity will not be local in the conventional sense [11]. Such effects could be modeled 
by a non-commutative space-time. 

Also in condensed matter physics the concept of non-commutative space-time 
is applied. The theory of electrons in a magnetic field projected to the lowest Lan- 
dau level can be naturally described by a non-commutative field theory [12-16], 
where |6/ijy| is again proportional to the inverse magnetic field. These ideas are 
relevant for the quantum Hall effect. For the integer quantum Hall effect, a non- 
commutative treatment serves as an alternative description to standard condensed 
matter techniques. This is already remarkable, since it is the first application of non- 
commutativity geometry that provides phenomenological results. However, with 
these methods only the known results are reproduced; it does not provide new in- 
sight in the nature of the integer quantum Hall effect. 

This may be different in the case of the fractional quantum Hall effect. That effect 
is not well understood from the theoretical point of view. Here a non-commutative 
field theory is considered by many researchers as the most promising candidate for 
a description. 

One may also try to study the non-commutative analog of pure Yang-Mills the- 
ories or of QED and QCD. Such theories can be considered as an extended standard 
model and a study of them could allow for an experimental verification or falsifi- 
cation. Hence results from these extended theories may provide sensible tests of a 
quantized space-time. 

The above described applications of non-commutative field theory suffer from 
the ultraviolet /infrared (UV/IR) mixing. This effect causes still severe problems in 
a perturbative treatment. Our goal is to study non-commutative field theories on a 
completely non-perturbative level. 

This work represents the first non-perturbative study of non-commutative field 
theories. As usual when entering a new topic we studied toy models, which share 
important properties of the full theory. In this work we studied field theory in 
lower dimensions and we focused on basic properties of these theories. The results 
presented here are published in Refs. [17-19]. 

In a two dimensional gauge theory we address the problem of renormalizabil- 
ity. It is an interesting question whether this model can be renormalized non- 
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perturbatively. In the three dimensional theory we studied the effects of UV/IR 
mixing. Our main interest was here the phase diagram of this theory and the ques- 
tion if there is a phase with spontaneously broken translation invariance, as it had 
been conjectured from analytic results. In addition we studied the dispersion rela- 
tion in this theory, for which perturbation theory predicts a deformation due to the 
non-commutativity. 

This work is organized as follows: in Chapter 2 we give an introduction to non- 
commutative field theory. We concentrate here on the main differences compared 
to commutative field theory. In a first Section we set up the non-commutative 
geometry on which we define a scalar field theory and a pure gauge theory (Sections 
2.2 and 2.3, respectively). These are the theories studied in this work. In Section 
2.4 we briefly comment on the extension to the non-commutative standard model. 
In addition to the motivation already given in the introduction we want to motivate 
the study of non-commutative space-time from this point of view. 

Chapter 3 is dedicated to the lattice regularization. As we already mentioned in 
this introduction the non-commutativity of space-time does not cure the ultraviolet 
divergences, and therefore one still has to regularize the theory. To this end we will 
introduce a momentum cut-off via discretization of space-time. 

In Chapter 4 and 5 we present the two models we investigated; a two dimensional 
non-commutative pure gauge theory and a three dimensional scalar field theory. The 
explicit construction of the lattice actions as well as the Monte Carlo results of our 
studies are presented in these Chapters. In Chapter 6 we show results on a two 
dimensional scalar theory, and in Chapter 7 we summarize our results and give an 
outlook. For the sake of continuity the details of the simulations are banned to 
Appendix A. 

Note that throughout this work we always work in Euclidean space-time. We 
should mention here that in contrast to the commutative case, it is an open question 
if the Euclidean version of non-commutative field theory can be interpreted in the 
Minkowski world, since there is no equivalent of the Osterwalder-Schrader theorem 
[20]. However, in non-commutative field theory with a commuting time coordinate it 
is generally beheved that this interpretation exists. For a discussion of Wick rotation 
and the related question of unitarity, see e.g. Refs. [21-23]. 
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Chapter 2 

Non— commutative field theory 



In this Chapter we give an introdnction to the concept of non-commutative space- 
time and field theories defined on it. We will work out the main differences to their 
commutative counterparts and discuss the additional problems that arise in such 
theories. For a general discussion of non-commutative space-time see for example 
Refs. [11,24]. We follow here the discussion in Refs. [25,26]. 

2.1 Non— commutative flat space— time 

In this Section we discuss those features of non-commutative geometry, which are 
needed to define field theories on such a geometry. We will find two alternative 
formulations, in terms of Weyl operators and in terms of functions with a deformed 
multiplication. 

2.1.1 Weyl operators 

Let us start with the commutative algebra of complex valued functions on d di- 
mensional Euclidean space-time W^. An element of this algebra corresponds to a 
configuration of a scalar field, with pointwise addition and multiplication. We con- 
sider here functions of sufficiently rapid decrease at infinity, so that any function 
f{x) may be described by its Fourier transform 



A non-commutative space-time can be defined by replacing the local coordinates 
x e M*^ by Hermitian operators X/^ satisfying 



The non-commutativity tensor 0^i, is antisymmetric with the dimension length 
squared and it can in general depend on space-time. Here we restrict ourselves 
to a constant non-commutativity tensor parametrized by the non-commutativity 




(2.1) 



(2.2) 
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= J) «) . (2.3) 

Here we assume the space-time dimension d to be even. The i;^ generate a non- 
commutative and associative algebra of operators. Elements of this algebra, the 
Weyl operators can be constructed by a formal Fourier transform involving 

the operators and the ordinary Fourier transform of f{x) 

W[f] = / ^ m e^'^'^ . (2.4) 

Combining equations (2.1) and (2.4) we find an explicit map A{x) between operators 
and fields 

W[f]=f(fxf{x)A{x) with A{x) = f -^e'^'^^^e-''''^''^ , (2.5) 
J J (27r) 

where A (a;) is a Hermitian operator that can be understood as a mixed basis for 
operators of fields. 

We may define a hnear and anti-Hermitian derivative on the algebra of Weyl 
operators by the commutator relations 

[^M, x^] = 5^y, [4, 4] = ic^y , (2.6) 

where c^^ is a real valued c-number. With this definition of the derivative one can 
show that 

[d,,^{x)] = -d,^{x). (2.7) 

Together with equation (2.5) and integration by parts one obtains that the deriva- 
tive of Weyl operators is equal to the Weyl operator of the usual derivative of the 
functions 

[4, W[f]] = J d'x d,f{x)A{x) = W[d,f] . (2.8) 
Any global translation x + v with v e R*^ can be obtained with the unitary operators 

A{x + v)^ e''^^^A{x)e-''-'^-' . (2.9) 

This follows directly from the commutator relation (2.7). This property implies that 

any trace TrA(a;), with Tr defined on the algebra of Weyl operators, is independent 
of a; G M°'. Together with equation (2.7) it follows that the trace Tr is unambiguously 
given by an integration over space-time 

TrW[f] = J d'^xfix), (2.10) 

with the normalization TrA(x) = 1. 
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In Ref. [26] it is shown that if is invertible (which imphcs that the dimension 
d of space-time has to be even) the product of two operators A at distinct points 
can be computed as follows 

A{x)A{y) = ^ I ^ I d'^zA{z)e-'''^^~'^>^"^''-'^^^y-'^" . (2.11) 
TT** det B J 

Together with the normalization of the trace the operators A{x) form an orthonormal 

set for X e M'', 

Ti{A{x)A{y))=5\x-y). (2.12) 
With this property of A{x) we can define the inverse map to (2.5) 

/(a;) = Tr(l^[/]A(x)). (2.13) 

This one-to-one correspondence can be thought of as an analog of the operator-state 
correspondence in quantum mechanics. 



2.1.2 The star— product 

Let us now consider the product of two Weyl operators corresponding 
to the two functions f{xi) and g{x2)- We want to transform this product into the 
coordinate space representation with the help of the inverse map (2.13), 

h{x)^l^{W[f]W[g\A{x)) . (2.14) 

To achieve this we rewrite the product in terms of the map (2.5) 

W'l/]»'l9l = y d'xi j (i''x2/(xi)9(x2)A(xi)A(x2) 



X f{xi)g{x2)A{x^)e 



where in the second line equation (2.11) was used. Multiplying both sides with A{x) 
from the right and taking the trace leads to 



h{x) ^TT{W[f]W[g]A{x)) 

^ / ^'"^ / ^'^2/(xO^(x2)e-^(«-^)-(--^)''(--^)^ ^2.16) 

= f{x) exp Qa^e^^S^ g{x) = f{x) -k g{x) , 

where we used the completeness relation (2.12). Using this product we obtain 

W[f]W[g] = W[f^g], (2.17) 
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i.e. the product of Weyl operators is equal to the Weyl operator of the star-products 
of functions in coordinate space. With the star-product we can rewrite the commu- 
tation relation between space-time operators (2.2) in terms space-time coordinates 

\p^fn •^I'li, •^^J. ^ •^i' ■^t' ^ "^A* ^0/1^ • (2. 18) 

The star-product is associative but non-commutative. For 0^j, = it reduces 
to the ordinary product of functions. It can be thought of as a deformation of the 
algebra of functions on R'^ to a non-commutative algebra, with the same elements 
and addition law, but with a different multiplication law given by (2.16). Note that 
the commutator of a function f{x) with the coordinates can be used to generate 
derivatives 

-k f{x) - f{x) 'kXf, = iQf,uduf{x) . (2.19) 
Due to the cyclicity of the trace defined in equation (2.10) the integral 

Tr{W[h]W[h]...W[U])^ J d''xMx)*Mx)i.---i.Ux) (2.20) 

is invariant under cyclic permutations of the functions fi{x). In particular the inte- 
gral of the star-product of two functions is identical to the integral of the ordinary 
product of two functions 

J dVi(^) */2(^) = J d''xh{x)h{x) . (2.21) 
We have now two ways to encode non-commutative space-time: 

• we can use ordinary products in the non-commutative C*-algebra of Weyl 
operators, 

• or we may deform the ordinary product of the commutative C*-algebra of 
functions in to the non-commutative star-product. 



2.1.3 The non— commutative torus 

In this Subsection we briefly discuss the case when space-time is a d-dimensional 
torus T'^ instead of W^. For a more detailed discussion, see [25,26]. 
Let us consider functions f{x) on a periodic torus 

f{x + T.^,fi) = f{x) with z/=l,...,d. (2.22) 

(x is the unit vector in the direction and S^,y is the d x d period matrix of the 
torus. Due to this periodicity the momenta kf^ in (2.1) are discretized according to 

A;^ = 27r (E-^)^^m^ with m^eZ. (2.23) 

Using the discrete version of Fourier transform (2.1) wc can define a mapping from 
fields to operators in the same way as we did in flat M*^. The result is 

A{x) = ^ ^ I JJ (^Z^^'^j ^-niE.<,^''pm.mp^-2niii::-^),^m.xi, ^ ^2.24) 



ki/=i 
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where we introduced the dimensionless non-commutativity tensor O 

e,. = 27r(E-^)^^e,.(E-^)^^ (2.25) 

and the operators 

^ g2^i(s-i).^x^ with Z^Z^ = e-2''*®'"'Z^Z^. (2.26) 

With this mapping we find the one-to-one correspondence (2.5) and (2.13) also on 
the non-commutative torus. These definitions will reappear in Section 3.1, where 
we discuss the discrete torus. 

2.2 Non— commutative scalar field theory 

Having defined the algebra of functions of non-commutative space-time we are able 
to define a scalar field theory on this geometry. At this point we make a change of 
notation and introduce the short-hand notation for Weyl operators — > /. 

2.2.1 Non— commutative scalar action 

We start our discussion with the action of an Euclidean commutative theory, 

S{^\ = y d'x Q d,<\>{x)d,<\>{x) + ^02(x) + \<\>\x)^ , (2.27) 
where is a real valued scalar field and d is the dimension of space-time. 

To transform an ordinary scalar field theory to a non-commutative field theory 
we can use one of the procedures described in the last Section. Either we may use 
the Weyl quantization via Hermitian operators 0, or we use the deformation of the 
product into the star-product (2.16). 

The quantum field theory written in terms of Weyl operators 0, corresponding 
to a real scalar field 0(a;) on M"*, reads 

I #exp(-5[0]) 

/I 2 X X (2.28) 

^[0] = Tr(^i[4,0]^ + ^0^ + ^0^J , 

where the kinetic term is a direct consequence of equation (2.8) (it involves a sum 
over \x). The measure is here the ordinary path integral measure for scalar fields 

This theory may be formulated in coordinate space by applying the map (2.13) 
to the action (2.28) and using equation (2.17), 



m = / 



d'^x 



^ {d^<^{xyf + "^(^{xf + ^0(x) * <^{x) * 0(x) * 0(x) 



(2.29) 
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the kinetic term and the mass term do not contain the star-product, because of the 
property (2.21) of the star-product. As a consequence the commutative and the 
non-commutative theory coincide for free fields. The difference arises due to the 
self-interaction term 




Tr(0^) = / d'^x (j){x) ★ (j){x) ★ (p{x) ★ (j){x) 



^ \ \ \ (2-30) 



.a=l / \a=l 



with the interaction vertex V in momentum space 

V{ki,k2,k3,h) ^ e-i^^'-^-^O'''^^''' . (2.31) 

This vertex contains a momentum dependent phase factor and it is therefore non- 
local. 

2.2.2 UV/IR mixing 

We discuss this important difference of non-commutative field theories compared to 
commutative theories at the example of one loop mass renormalization of the 4d 
theory, given by equation (2.29). To this end we consider the one particle irreducible 
two-point function 

oo 

r(p) = {4>{p)4>{-p)) = E A"r(")(p) ■ (2-32) 

n=0 

At lowest order the two-point function is given by r(o)(p) = + m? . The one 
loop contribution splits topologically into two parts, one planar and one non-planar 
diagram 

= 3/ (27r)4P + m2' ^^-^^^ 



r(i)^^^ - ^ / ^ exp 



1 f d'^k exp {ik^p^Q ^y) 
6 J (27r)4 P + m2 

In Refs. [3,4] it is shown that the contribution of planar diagrams to non-commutative 
perturbation theory is proportional to the commutative case (to all orders). There- 
fore the planar divergences may be absorbed into the bare parameters, if and only if 
the corresponding commutative theory is renormalizable. This already disproves the 
expectation that non-commutative quantum field theory would not require renor- 
malization. 

In the case of the non-planar diagrams the situation is different. Rewriting the 
denominator in equation (2.34) in terms of a Schwinger parameter 

-— ^= / dae-^(''+^'\ (2.35) 
k^ + Jo 
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planar non-planar 

Figure 2.1: The planar and non-planar one loop contribution to two-point function 
(2.32). 

and introducing a momentum cut-off A by multiplying the resulting integrand in 
(2.34) with a Pauli-Villars regulator exp(— 1/(A^q;)), leads in d dimensions to [26] 



where Kn{x) is the irregular modified Bessel function of order n. In d = 4 the leading 
divergences of equation (2.36) are given by 

= ^ (Ae\ - log (^) + 0(1)) . (2.37) 

Here we introduced the effective cut-off Aeff given by 

^eflf = 7^ = 1 , 02 2 ■ (2.38) 

The two-point function rip remains finite in the limit A — > 00, because it is ef- 
fectively regulated by the non-commutative space-time. The complete one loop 
corrected propagator then reads 

r(p) = p2 ^ + 2Ar«(o) + ArW(p) + o{x') , (2.39) 

where we used 

r?^ = 2rS(p = o). (2.40) 

The UV limit (A oo) does not commute with the IR limit (p — > 0) or with the limit 
6^0. At small momenta or small non-commutativity parameter the two-point 
function reads 

Tip) ~ p2 ^ + 3\T[^iO) + 0(A2) . (2.41) 

Taking now the UV limit leads to the standard mass renormalization of the 
theory. Taking these limits vice versa, the effective cut-off is given by 

Ai=^ (2.42) 

and Aeff diverges — and therefore also Tnp{p) — either in the limit ^ ^ — > or in 
the infrared limit when the incoming momentum p goes to zero. We may absorb 

^Note that the limit ^ ^ in the non-commutative action (2.29) leads to the commutative 

action, since in this limit the star -product turns into the usual product. In the quantized theory, 
after the cut-off A is removed, the limit 9^0 does not lead to the commutative theory, as equation 
(2.43) shows. 
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the planar one loop contribution of T{p) by defining the renormalized mass through 
= + 2Arnp (0). Removing the cut-ofi^ while keeping M^^ fixed, then leads to 
a finite r(p) for finite incoming momenta p. For zero momentum r(p) diverges and 
the divergence at one loop is given by 



with ^ = gg^- Here we see that a non-zero non-commutativity tensor O^j, replaces 
the standard ultraviolet divergence with a singular infrared behavior. This mixing 
between high and low energy effects is called UV/IR mixing. 

The long distance behavior of the spatial correlators is controlled by the pole in 
the upper half plane nearest to the real axis, as in the commutative case [11]. Due 
to the additional term in equation (2.43) the poles of the propagator are now at 



In the weak coupling limit this pole is here at p = i\fU^/ (M^e 6*) and not at p = iM^^ 
as in the commutative case. This can be interpreted a a new mode with mass 
1712 = v^^/(Meff^). By definition m2 is much smaller than Afcff and therefore 
dominates the behavior of the spatial correlators. In the commutative theory these 
correlation functions decay exponentially if Meff > 0. Here we obtain at small A a 
power-like decay of the correlators, leading to the long range correlations [4]. At 
large A the decay is again exponential, but now with a decay constant oc -\/A. 

2.2.3 Phase structure of non— commutative theory 

The UV/IR mixing is one of the most interesting properties of non-commutative 
field theory and has no counterpart in the commutative case. A number of new 
effects and also problems occur due to this term. In particular the phase diagram of 
the non-commutative X(f)^ model is changed, which we want to discuss here. 

The low momentum singularity of r(p), discussed in the last Subsection, has 
a large impact on the phase diagram. Since a phase transition should involve the 

momenta that minimize T{p), it is not likely that the low momentum modes will 
participate in a phase transition. If there is a phase transition at all it will be 
driven by non-zero momentum modes. Then the IR divergence leads to an oscilla- 
tion in the sign of the correlator (0(0) ★0(a;)), indicating a new type of ordered phase. 

Gubser and Sondhi studied the phase diagram of 4d A0'* theory [27] within the 
framework of a one loop self-consistent Hartree-Fock approximation [28]. We do 
not discuss their calculation, but summarize their results on the phase diagram. 

At small non-commutativity parameter 9 they obtained an Ising type (second 
order) phase transition leading to an uniformly ordered phase with {(f)) ^ 0. At 



r(p)=p^ + M,^ff + e 



A 



+ subleading terms 



(2.43) 




(2.44) 
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Figure 2.2: The phase diagram of non-commutative 4d Xcf)'^ theory conjectured by 
Gubser and Sondhi [27] in the m? - X plane. 



sufficiently large 9 the minimum of r(p) is not at p = 0. The phase transition is 
now driven by modes p 7^ and it is of ffist order. This leads to an ordered phase 
where the translation invariance is broken spontaneously. In this phase ((/>) varies 
in space, which implies the ground state to involve some non-uniform patterns like 
stripes. These patterns depend on the momentum mode which drives the phase 
transition. In Ref. [27] these results are summarized in a qualitative phase diagram 
in the ra^hr^-OK^ plane, where A is a momentum cut-off. 

In Section 2.2.1 we have seen that only the interaction term depends on the 
parameter 9. Therefore increasing the coupling A also amplifies effects of non- 
commutativity. According to Ref. [27] the phase diagram in the m?K~^-\ plane is 
then given by Figure 2.2. 

A similar phase structure was conjectured in three dimensions. In two dimensions 
it was argued that a striped phase does not occur. Gubser and Sondhi worked with 
an action of the Brazovskiian form [28], which is local. Hence the Mermin- Wagner 
theorem [29-31] applies, which states that spontaneous breakdown of a continuous 
symmetry is not possible in two dimensions. We come back to this point in Chapter 
6. 



In another approach renormalization group techniques were used to study the 
phase diagram of the model [32]. Chen and Wu obtained in d = 4 — e a new 
IR stable fixed point, i.e. the non-commutative counterpart of the Wilson-Fisher 
fixed point. This fixed point is stable, and therefore a striped phase exists, when 
9 > 12/y^. In contrast to the results in Ref. [27], this implies that in li = 4 there is 
no striped phase. Since we studied the 3d model, we will not address this controversy 
in this work. 
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2.3 Non— commutative gauge theory 

In this Section we extend our considerations to gauge theories defined on a non- 
commutative plane. 



2.3.1 Star— gauge invariant action 

To define a Yang-Mills theory on a non-commutative plane we have to generalize 
the map (2.5) in Section 2.1 to the algebra of n x n matrix valued functions. Let 
A^{x) be a Hermitian gauge field on W^, which corresponds to the unitary gauge 
group U (n). We can introduce the Weyl operators corresponding to A^i^x) by taking 
the trace of the tensor product of A(a;) and the gauge field 



(2.45) 



where A{x) is defined in equation (2.5). Based on equation (2.8) a non-commutative 
version of the Yang-Mills action can be defined 



S[A] 



1 



-TrtrTv 



dfi, Ajy 



du, A^ — i A^, Ay 1 , 



(2.46) 



where the term in brackets is the operator analog of the field strength tensor. Here 
Tr is the operator trace (2.10) and tr^r denotes the trace in color space. This action 
is invariant under transformations of the form 



A^^GA^G^-iG 



(2.47) 



where G is an arbitrary unitary element of the algebra of matrix valued operators, 



1. e. 



GG^ = G^G = i®ln. 



:2.48l 



The symbol 1 is here the identity on the ordinary Weyl algebra and is a n x n 
unit matrix. 



To set up the action in coordinate space we can construct an inverse map of (2.45). 
By mapping the product of matrix valued Weyl operators to coordinate space, using 
this inverse map, again the star-product (2.16) appears. The Yang-Mills action in 
coordinate space then reads 

S[A] = J d'^xtrjv (F^,(x) *F^,(x)) , (2.49) 

where we introduced the non-commutative field strength tensor given by 

F^y{x) = d^Ay{x) - dyA^{x) - i [A^{x),Ay{x)l . (2.50) 

The index indicates that the products in this commutator are star-products. 
Prom equation (2.50) we see that already for the simple gauge group U{1) we have 
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a Yang-Mills type structure. Therefore there exist three and four point gauge inter- 
actions and non-commutative U{1) theory is asymptotically free. 

The invariance under unitary transformations in operator space translates here 
into an invariance of the action (2.49) under star-gauge transformations given by 

A^{x) G{x) ★ Af,{x) ^ G{xy - iG{x) ^ d^G{xy , (2.51) 

where G{x) is a star-unitary matrix field, 

G{x)i.G{x)'' = G{x)UG{x) = In. (2.52) 

Equation (2.52) is equivalent to the unitarity condition (2.48). 

So far we considered non-commutative U{n) theories which reduce to the or- 
dinary U{n) theories in the limit ^ — > (on the classical level). In Ref. [33] it 
was shown that for other gauge groups like SU{n) this cannot be realized on non- 
commutative flat space. ^ The U{n) group is closed under the star-product; the 
product of two star-unitary matrix fields is again star-unitary. In contrast to U (n) 
the special unitary group SU (n) is not closed, since in general 

det(G)*det(ii') ^det(G'*ii'). (2.53) 

The U{1) and the SU{n) sectors in the decomposition 

U{n) = U{1) ® SU{n) /Z„ (2.54) 

do not decouple in the non-commutative case, because the U{1) photon interacts 
with the SU{n) gluons [34]. 

2.3.2 Star— gauge invariant observables 

To construct star-gauge invariant observables we consider an arbitrary oriented 
smooth contour Cy in space-time, which connects the points x and x -\- v. The 
holonomy of the gauge field over this contour is described by the non-commutative 
parallel transporter 

U{x- C„) = P exp^ J di^A^ix + 0^ , (2.55) 

where P indicates path ordering and ^ parameterizes the contour. The index at 
the exponential function indicates that in the expansion of this function the star- 
product has to be used. The parallel transporter (2.55) is a n x n star-unitary matrix 
field and transforms under the star-gauge transformation (2.51) like 

U(x;Cy) G(x)i.U(x;Cy)i.G\x + v) . (2.56) 

^We refer to a constant non-commutativity tensor 
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A remarkable fact in non-commutative field theory is that space-time translations 
can be arranged by (star-) multiplication with plane waves 



G{x + v)^ e^'^''^" * G{x) ★ e-''""'''' with A;^ = (e^^) ^ , (2.57) 



where we assume to be invertible. That this equation holds can easily be shown by 
expanding the exponential functions and using equation (2.19). With the definition 
of the non-commutative parallel transporter and equation (2.57) we can associate a 
star-gauge invariant observable with any arbitrary contour Cy by 



It is straightforward to show the invariance under the star-gauge transformation 
(2.56) by using equation (2.57) and the cyclicity of the trace over the star-product. 

In commutative gauge theory gauge invariant observables can only be constructed 
from closed loops. In contrast to that, equation (2.58) shows that in non-commuta- 
tive gauge theory we can find star-gauge invariant observables associated with open 
contours. The vector k in equation (2.58) can be regarded as the total momentum 
of the open loop. This is again a manifestation of the UV/IR mixing phenomenon, 
discussed in Section 2.2.2. If we increase the momentum in a given direction, the 
contour will extent in the other directions according to Q^^ku- 

This completes our introduction to non-commutative field theories in the con- 
tinuum. We showed how to define a scalar field theory and a pure gauge theory 
on a non-commutative plane, and we discussed the main differences compared to 
the commutative case. This sets up the framework for our numerical studies to be 
presented in Chapters 4 and 5. In the next Section we will discuss further properties 
and problems of non-commutative field theory. 



2.4 Phenomenological implications of a quantized 
space— time 



In this Section we present phenomenological consequences of a ^-deformed space- 
time. To this end we discuss briefiy some aspects of the non-commutative standard 
model. 

Gauge fields coupled to matter fields 

To set up the non-commutative standard model we have to extend our considerations 
in Section 2.3 to the case where the gauge field couples to matter fields. We start 
from the action of a free Dirac field 




(2.58) 




(2.59) 
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where we extended the commutative theory to a non-commutative theory by re- 
placing the usual products of fields with star-products. The Grassmann valued 
fermionic fields are represented by ^(a;). To obtain an action that is invariant under 
the star-gauge transformations 

il){x) ^ G{x)i.ij{x) and ^(x) ^ ^(x) ^ G^(a;) , (2.60) 

with a star-unitary nxn matrix field G{x) , we have to modify the kinetic term of the 
action (2.59). We follow here Ref. [35] and introduce, in analogy to the commutative 
case, the covariant derivative 

Df,ilj{x) = d^i}{x) + igA^{x) ★ V'(x) , (2.61) 

where A^{x) is the gauge field that generates the unitary group U (n). The derivative 
(2.61) is covariant under the star-gauge transformation (2.51), 

D^^G{x)i<D^i<G\x). (2.62) 

Replacing the derivative in (2.59) by the covariant derivative leads to the star-gauge 
invariant fermion action 

^fermion[^, i'M^ j d'^^ {lp ★ [Tm^M + m] ^j) . (2.63) 

The complete action is then given by the sum of the gauge action (2.49) and the 
fermion action (2.63) 

S[A, ^, ^] = Syu[A] + 5fe™[^, ^, ^] ■ (2.64) 

As we mentioned in the last Section, it was so far not possible to formulate a 
non-commutative field theory for the special unitary group SU (n). The gauge group 
is restricted to U{n), which is (for n > 1) not a gauge group of the standard model. 
There are attempts to modify the space-time non-commutativity in order to take 
also SU{n) gauge groups into account [36,37], but this is an ongoing field of research. 
However, we may consider the star-unitary U{1) gauge field coupled to fermions as 
an extension of commutative QED. 

About renormalizability 

The UV/IR mixing poses severe problems for the renormalization of perturbation 
theory, which are not overcome yet. For instance, scalar fields can become unstable 
— tachionic — due to IR effects as the non-commutativity is switched on. The 
attempts to renormalize perturbation theory include methods known from standard 
field theory [38], the application of Wilson's renormalization group technique [39-41], 
controlling the IR divergences in the framework of supersymmetry [42,43] and the 
Hartree-Fock method [27]. 

In spite of some plausibility arguments in favor of perturbative renormalizability, 
it is an open question if non-commutative quantum field theories do really have 
finite UV and IR limits. There are even claims that basic non-commutative field 
theories, like non-commutative QED, are not renormalizable [44]. 
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Violation of Lorentz symmetry 

Since 0^,^ carries Lorentz indices, two distinct types of Lorentz symmetries have to 
be considered [45]: tlie observers Lorentz transformation and tlie particle Lorentz 
transformation. Tlie action (2.64) as well as the scalar action (2.29) are fully covari- 
ant under rotations or boosts of the observers reference frame, because 0^1, and the 
fields transform covariantly. This does not hold anymore in the case of rotations or 
boosts of a particle [46]; ©^^ is unaffected by these transformations. ^ 

The broken Lorentz symmetry implies a deformed dispersion relation, where 
the deformation depends on the theory under consideration. In the case of the non- 
commutative theory the dispersion relation is given by the poles of the irreducible 
two-point function (2.43). We focus here on the case with two commuting and two 
non-commuting directions. The on-shell condition then reads 

E{p)' ~ + P' + Ml^ + , (2.65) 

where p is the momentum that corresponds to the non-commutative plane and P 
corresponds to the commutative space coordinate. We include here only the leading 
IR divergence in first order of A; in addition there is also a logarithmic divergence. 

Phenomenology of non commutative space— time 

The effects of non-commutativity on quantities that are measurable in experiment 
are studied intensively by many authors. These studies may be interpreted as an 
attempt to set stringent limits on the non-commutativity parameter 0, or to suggest 
possibly measurable effects of non-commutative space-time. We will present here 
an incomplete list of current investigation in this field. 

Deformed dispersion relation in gauge theories 

Also in gauge theories Lorentz symmetry is broken, which leads to a deformed 
dispersion relation and changes the particle propagation. This was studied for exam- 
ple in Refs. [47,48]. There a dispersion relation for photons was obtained containing 
a term as in the scalar case. Based on one loop perturbation theory the authors 
of Ref. [47] obtained the photon dispersion relation 

E{p)^ = p^ + , (2.66) 

where they chose the time to be commutative. The coefficient depends on the 
number of bosonic and fermionic degrees of freedom present in the theory and g is 
the coupling constant. 

This may give rise to experimental verification of quantized space-time. The 
non-linear term in E{pY leads to a vacuum dispersion of light, such that the speed 

^Thc non-commutativity tensor O^^ plays the role of an inverse background field (see introduc- 
tion) . From this point of view the Lorentz symmetry breaking occurs very naturally. The symmetry 
would be restored in a formalism that transforms the particle fields as well as the background field. 



2.4 Phenomenological implications of a quantized space— time 21 



of light depends on the wave length. In addition the dispersion depends on the 
polarization of the photons leading to birefringence. 

The authors of Ref. [49] suggest to measure this deformation of the dispersion 
relation by so-called time of arrival measurements. In these measurements the time 
delay between photons with different wave length, emitted simultaneously up to a 
known precision, is measured. The delay will depend on the time the photons are 
traveling. This effect is (if existent at all) very small. Therefore the sources of the 
photons should be far distant to accumulate the delay to a measurable effect. 

There are attempts by experimentahsts to study this effect. Already in 1998 high 
precision measurements ^ of an energy dependent speed of light were performed [50] . 
However, up to a given sensitivity these measurements did not show any energy 
dependence. 

Newer experiments might give more insight. For example the HESS project [51], 
just started to take data, has a higher sensitivity than in Ref. [50]. Besides other 
projects, this collaboration intend to make measurements related to the vacuum dis- 
persion [52]. Another candidate for measuring such effects is GLAST [53], which is 
expected to start measurements by the year 2005. 



Threshold anomalies 



It is still an open question why high energy photons {E > 20TeV) from far distant 
galaxies can be detected on earth [54]. Photons in this energy range traveling over 
galactic distances should interact with the cosmic microwave background, producing 
electron positron pairs 

77 — > e"^e~ . 

The threshold for this reaction is approximately E pa ITeV [55, 56] and it should 
make the observation of high energy photons very unlikely. 

A second threshold exist for cosmic high energy protons. There the interaction 
of the protons with the cosmic microwave background leads to 

p + 7 p + n . 

The protons loose energy when producing pions and the threshold, i.e. the GZK 
threshold, for this process is at approximately E ^ 5 x 10^^ eV [55,57] Therefore 
cosmic protons with higher energy should not be observed. However, cosmic proton 
rays have been detected beyond this limit 



In both cases non-commutative space-time could provide an explanation of these 
anomalies. The thresholds can be estimated kinematically using the momentum and 
energy conservation. In a Lorcntz invariant theory the threshold momentum in the 
case of high energy protons is approximately [59] 

nipm^ + ml 

Pthreshold ~ , (2.67) 



*The accuracy in this measurements was 6T/T « 10 where ST is the time delay and T is 
the overall time. 
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where rup and m^r are the masses of the proton and of the pion, respectively, and 
is the energy of the background photon that interacts with the proton. In a 
non-commutative space-time the dispersion relation is given by equation (2.66). In- 
serting the deformed dispersion relation into this computation leads to a ^ dependent 
threshold. 

Bounds on 9 

So far we discussed new effects and known problems that might be explained by 
non-vanishing 6. Let us now address the limits on 6 set by existing experiments and 
the possibility of setting bounds on 6 in near future experiments. 

The fact that in the aforementioned time of arrival measurements no energy de- 
pendence was measured, allows to set a lower limit on 9 (if we assume a non-zero 
9 from the beginning). Note that here we find a lower limit, because the dispersion 
relation (2.66) implies that larger values of 9 correspond to a softer deformation of 
the standard dispersion. ^ On the other hand if one wants to solve the cosmic proton 
threshold anomaly within the framework of non-commutative field theory one finds 
an upper bound of 9. According to Ref. [47] this leads to a rough estimation of the 
a parameter 9 in the range (lO^TeV)-^ <9 < (lOTeV)-^ 

More precise bounds on 9 could be obtained from accelerator physics. Since non- 
commutative QCD is not well formulated yet, these investigations are restricted to 
measurements described by QED. Hence one expects the most significant results 
from linear colliders. 

The idea is to compute for example cross sections for scattering processes from 
the standard model and from the non-commutative standard model [60]. ® This 
leads to different predictions for the two models, where both predictions depend of 
course on the center of mass energy ^/s. Measurements of these observables then 
may verify the non-commutative picture or set an upper bound on 9. Measure- 
ments in existing colliders did not show any signal of non-commutative space-time. 
Therefore one is looking for search limits, which can be probed in future colliders 
(TESLA [64]. NLC [65] and CLIC [66]) 9. For example. Ref. [67] reports search 



process 


search limit on l/s/O at -/s = 0.5 — 5 TeV 


pair production 7 7 — > e~ e+ 


500 - 2700 GeV 


Compton scattering 67^67 


1000 - 6500 GeV 



Table 2.1: Sensitivity to measure effects of non-commutativity in future colliders. 



limits for the processes hsted in Table 2.4 on the basis of the design of these coUid- 

■''Thc rather unexpected lower bound of 6 is related to the order of UV limit and commutative 
limit, as we discussed in Section 2.2.2. This lower bound can be understood in the sense that if 
there is a non-commutative space-time, then 6 has to be larger than a certain value. However, 
this does not exclude 9 = 0. 

^The calculations in non-commutative QED rely on the Feynman rules developed in Ref. 
[34,61-63]. 
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ers. Further search hmits can be found for example in Ref. [68]. For recent reviews 
concerning the activities in this field see for example Ref. [69-71]. 

Note that the theories, which the above presented results are based on, still suffer 
from IR divergences. They are obtained from one loop calculations. Higher orders 
in perturbation theory are not under control yet. It is still an open question if the 
non-commutative field theories are renormalizable in two or three loop calculations. 
The limits on 9 were obtained assuming that higher order correction will not change 
the results qualitatively. However, one cannot exclude dramatic changes coming 
from higher loop contributions. 

2.5 Summary 

Let us briefly summarize the effects of non-commutative space-time on field the- 
ories defined on it. We showed that introducing non-commutativity via the com- 
mutation relation (2.2) leads to a non-commutative and non-local product, i.e. the 
star-product. Field theories may be formulated on this geometry by replacing all 
products in the commutative theory by the star-product of the fields. This results 
in a non-local action. In a perturbative treatment one discovers a mixing of high 
and low energy effects. The UV/IR mixing causes serious problems in perturbative 
renormalization, since some of the UV divergences in commutative theories turn 
into IR divergences in the non-commutative model even in the massive case. We 
discussed this new effect in one loop perturbation theory. Already there the diver- 
gences at low momenta cause enormous problems. The difficulties increase when 
these loops are sub-loops of a higher order contributions. 

In a ^-deformed space-time the Lorentz symmetry is explicitly broken. This 
leads to unexpected particle propagation like a momentum dependent speed of light. 
These effects are the most likely candidates for experimental verifications. 

In non-commutative gauge theories the gauge symmetry turns into a star-gauge 
symmetry. This allows us to construct star-gauge invariant observables that are 
associated with open Wilson loops. The open loops carry a momentum proportional 
to the separation between the endpoints. This is again a UV/IR mixing effect. In 
non-commutative QED the photons interact, which might again give rise for new 
measurable effects. 

As an alternative to perturbation theory we are studying non-commutative field 
theories in the lattice approach [72]. This work has to be considered as the beginning 
of a long term project. At the end of this project we want of course to compute 
phenomenological quantities, but in the starting phase we will study field theory in 
lower dimensions in order to understand better the non-perturbative treatment. 

However, already the study of these toy models gives insights into the four di- 
mensional theory. Since also for these models there exist perturbative results, it is 
an interesting question if in a non-perturbative study new effects arise and how far 
perturbation theory can be confirmed. 
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Chapter 3 

Lattice regularization 



As we have seen in Section 2.2.2, introducing non-commutativity of space-time does 
not cure the uhraviolet divergences occurring in the commutative case. Therefore 
the field theories still have to be regularized, by introducing a momentum cut-off. 
In this Chapter we discuss the lattice regularization. For a more detailed discussion 
we refer to Refs. [25,26,73-75]. 



3.1 Discrete non— commutative space— time 

On the cubic lattice the space-time points are restricted to discrete values G 
aZ, where a is the lattice spacing. Here momentum space is compact and the 
momenta A;^ must be identified under the shift 

27r 

ki^^kf^-\ with i/ = l,2,...,d. (3.1) 

a 

As in the continuum, non-commutativity comes into the game by replacing the 
ordinary coordinates by Hermitian coordinate operators X/^, which satisfy the 
commutator relation (2.2). As a consequence of equation (3.1), we can set up the 
operator identity 

By multiplying both sides of equation (3.2) with exp(— i/cpXp) we find 

—Xfj.Sfj,p j exp I Q,j,uki,—Si^p 1=1 with p=l,2,...,d, (3.3) 

where i is again the identity of the algebra of Weyl operators. The usual constraint 
of lattice field theory that the discretization has to be compatible with the spectrum 
of the position operator leads to 



2-n-i 



e— ^^^ = 1 for /X = l,2,...,d. (3.4) 
Moreover we find an additional constraint for the momenta 

e^^A;^ e 2aZ with = 1, 2, . . . , d . (3.5) 
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Combining this result with the periodicity in momentum space (3.1) imphes that 
Qn^n/a^ is an integer for all /x, leading to a discrete momentum 

k/j, — 2n (Ti~^^ ^^rrit, with rrit, & Z; i' — l,...,d. (3-6) 

The periods E^,^ are integer multiples of the lattice spacing a and the periodicity 
(3.1) can be expressed by the integers rrii, via rrii, — > rrii, + ^^^i^- 

Due to the discrete momentum space the space-time coordinates are restricted 
to a periodic lattice 

+ E^^ with i/ = 1, . . . d . (3.7) 

Therefore the lattice regularization forces space-time to be compact. The discrete 
compactification is a non-perturbative manifestation of the UV/IR mixing effect; the 
perturbative manifestation was described in Section 2.2.2. For practical applications 
this means that systematic errors resulting from the finite lattice spacing are always 
entangled with finite volume errors. 

The continuum limit a — > does not commute with the commutative limit 
|0| 0. Taking first the continuum limit restores the infinite space-time IR'^, while 
taking first the commutative limit shrinks the space-time lattice to a single point as 
one infers from equation (3.5). 

Due to the momentum discretization (3.6) we cannot use the unbounded oper- 
ators Xfj,. Instead we have to use the unitary coordinate operators Z^, which have 
been introduced in Section 2.1.3 

= e'''^^^"').'/'' . (3.8) 

These operators generate the algebra of functions on a non-commutative torus ^ 

4z, = e-2-®-Z,Z^, (3.9) 

where the dimensionless non-commutativity tensor is defined as 

e,. = 27r(E-i)^^e,,(E-^)^^. (3.10) 

Moreover since we are dealing with a discrete torus, one should rather use the shift 
operator Z)^ 

i^^ = e"^-, (3.11) 

instead of the linear derivative (9^ (defined in (2.6)). This operator effects translations 
in units of the lattice spacing a 

D,Z^Dl^e''''<''~'^'^.Z^. (3.12) 
The algebra defined in (3.9) and (3.12) replaces the algebra based on (2.2) and (2.6). 



Note that in equation (3.9) and (3.12) no sum over n, v is taken. 
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Now we consider scalar functions f{x) on a non-commutative discrete torus, 
where we assume that these functions may be expressed by a discrete Fourier trans- 
formation of f{m). Lattice Weyl operators are then defined by the formal Fourier 
transform 

f^Y^ f{m) e2'^^(S"')''M"^.£M ^ (3.13) 

m 

analogous to the continuum (2.4). The relation between the integer vector m and 
the momentum k is given by equation (3.6). Replacing /(m) in equation (3.13) by 
its Fourier transform and using equation (3.9) leads to an explicit map A(a;) between 
operators and fields 

A(x) = 3^ ( JJ i^")] e~'''^-<p^'"''^"'^<'e~'^'''^^~^^''^'^"'''' . (3.14) 

This map is the lattice analog of the continuum map (2.24). Here x is a point on 
the space-time lattice (aZ)*^. The map between fields and operators then reads 

/ = ^/(x)A(x), (3.15) 

X 

fix) = Tr (/A(a:)) , (3.16) 
where the operator trace is uniquely given by 

Tr/ = J2 /(^) and TrA(x) - 1 for all x e (aZ)'^ . (3.17) 

X 

The lattice version of the star-product (2.16) can be found by applying the 
inverse map (3.16) to the product of two Weyl operators / g 

/(a;)*^(x)'='Tr(/^A(x)) ^ J]e-'<®"U''^7(^ + y^^^^ (3.18) 



y,z 



As in the continuum we obtain two equivalent descriptions of non-commutative 
space-time; either we use Weyl-opcrators or functions on (aZ)^ with a deformed 
product. 



3.2 Non— commutative field theory on the lattice 

With the definitions of the last Section it is straightforward to construct field theories 
on a non-commutative lattice. In the case of a scalar theory we use the map (3.15) 
to define the operators 4>. The action in terms of Weyl operators then reads 

= Tr (^^ E i^^^^l + + l^"^ > (3-19) 



28 



Chapter 3 Lattice reguleirization 



where -D^ is the lattice shift operator defined in equation (3.11) . Using the inverse 
map (3.16) leads to the lattice action of the field (j){x) 



X \ n 



X + afi)—(f){x)) H 4>'^{x) + —(f){x)-k((){x)-k(f){x)-k(l){x) 



2 ' ' ' 4 

(3.20) 

where fi is the unit vector in the direction /i. Since also the lattice version of the 
star-product is invariant under cyclic permutation, only the interaction term differs 
from the commutative case, as in the continuum. 

To construct a non-commutative Yang-Mills theory on the lattice we need the 
non-commutative counterpart of the link variables U^{x). In order to preserve star- 
gauge invariancc the link variables have to be star-unitary. This can be formally 
achieved by taking the parallel transporter defined in (2.55) between the lattice sites 
X and x + ap, 

(i-x+ajl \ 
I J d^A^iOj , (3.21) 

where is a non-compact U{n) gauge field and cxp^ indicates that in the expansion 
of the exponential function the star-product has to be used, t/^ is star-unitary and 
transforms under star-gauge transformations as 

U^{x) G{x) ★ U^{x) ★ G\x + ajl) , (3.22) 

where G{x) is again a star-unitary matrix field introduced in equation (2.52). The 
corresponding Weyl operators are given by 

C/^-5^A(x)®[/^(x). (3.23) 

X 

The star-unitarity of the link variables is now equivalent to unitarity of the Weyl 
operators. 

The action of a non-commutative Yang-Mills theory can again be constructed 
with the use of the lattice shift operators (3.11) or with the star-product 



S = J]Trtr^ [U, [d^U.DI) (dmIdI) Ul 



\ (3.24) 

where tr^v is the trace in the fundamental representation of the unitary group U (N). 

As in the continuum we find in coordinate space a simple recipe to construct 
a field theory on a non-commutative plane: replace all usual products in the com- 
mutative theory by star-products. In this spirit we now define star-gauge invariant 
observables in non-commutative lattice gauge theory, by replacing the products in 
the commutative Wilson lines U{x, C^) with star-products 

U{x, Cy) = U{x, III) ★ U{x + /ti, //2) * . . . . (3.25) 
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Here is an arbitrary lattice contour with the displacement vector v between the 
two end point of the contour. In the commutative case the trace over these lines 
is gauge invariant if the lines are closed, either on the plane or over the boundary. 
This implies v — or — E^j^n^, respectively, where e Z indicates how often 
the contour winds around the torus in the z/-th direction. 

On the non-commutative torus we can construct star-gauge invariant observables 
for any separation vector v by multiplication by plane waves 

0{C^) = ^ trjv C„) ★ 6^'='^^" , (3.26) 

X 

where the relation between the momentum k, carried by the contour, and the sepa- 
ration vector V on the torus is given by 

V,j. = Q^iuK + '^^nyfly . (3.27) 

The fact that open loops carry a momentum proportional to the separation vector 
V again displays the UV/IR mixing. 



3.3 Matrix model formulation 

In the last Section we constructed a scalar and a pure gauge action on a non- 
commutative lattice in terms of Weyl operators and in terms of functions with a 
deformed product, i.e. the star-product. The latter allows for an intuitive extension 
from commutative field theories to their non-commutative counterparts and it may 
seem to be suitable for a direct simulation. However, from a practical point of view 
this formulation rises enormous problems for Monte Carlo simulations. 

Already in the scalar case we have to implement the lattice version of the star- 
product, which involves a sum over the complete lattice for every product. This 
increases the needed computer time to an almost unreachable amount. In addition 
to this, in gauge theory we have to construct star-gauge invariant link variables U^, 
which makes the simulation even more expensive. To avoid these problems we will 
use a finite dimensional representation of the operator description. 

Here we are leaving the general discussion and restrict ourselves to the case that 
we studied numerically. This is a 2d non-commutative space or subspace with the 
period matrix given by 

E^, = NaS^, , (3.28) 

where a is the lattice spacing on a hyper-cubic lattice and N is the number of lattice 
sites in each direction. Then the coordinate operator (3.8) and the dimensionless 
non-commutativity tensor (3.10) are given by 

The quantization of momentum and the non-commutativity parameter 6 read 

= and O^l^'^N. (3.30) 
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The algebra (3.9) then simphfies to 

44 = e-^^'^'-^/^i.^ ; D^Z^bl = e^'^'^'^-'/^i, , (3.31) 

where e^i^ is the totally antisymmetric tensor. 

The key step is now to replace the operators and in the algebra (3.9) by 

N X N matrices, satisfying the relations (3.31). For odd values of N there is a simple 
choice to achieve this, given by the so-called twist eaters T^. In two dimensions this 
amounts to 

Di = ri = Si+ij modiV and D2 = T2 = Z^t^Si j 

= (r,)<-')/^ and z.^(rl) 

where we introduced the twist Z^^ in a way that will be useful in the next Section. 
Inserting (3.32) in (3.31) fixes the twist to 

Zu = ^21 = exp (7ri(A^ + 1)/N) . (3.33) 

The equations (3.32) are of course only one possible choice. For a general construc- 
tion see for example Refs. [25,73,76]. 

With these definitions an explicit form of the star-product (3.18) is provided and 
the operator A.{x) becomes a, N x N matrix 

N 

A(x) = ^mi^m2g-27rjmim2/iVg-2^im^x^/JV_ ^33^^ 

mi, 1712=1 

The map between the functions f{x) and the N x N matrices / in terms of A(x) 
then reads 

X 

where the trace now refers to the N x N matrices. If we insert A (a;) explicitly in the 
left equation of (3.35) we obtain the matrices J{m). Since the Fourier transform of 
f{x) is given by 

/(^) = ^ E /(^)e-''^^"^''^''/^ , (3.36) 

X 

the matrices J(m) read 

J(m) = z'P^z'^2^-'i^imim2/N ^337^ 
and we obtain a map between Weyl operators and the Fourier transform f{m). 

/ = ^ E fMAm) and /(m) = ^Tt (/ j\m)) . (3.38) 

m 

This completes the construction of a finite dimensional representation of field 
theory on a non-commutative lattice. In Chapter 4 and 5 we will construct explicitly 
the action of 2d Yang-Mills and 3d scalar field theory (with two non-commuting 
dimensions) in this representation. 



Chapter 4 

Numerical studies of 

non— commutative gauge theory 

In this Chapter we present the results of our studies of 2d non-commutative U{1) 
gauge theory. As aheady discussed in the previous Chapter the lattice version of 
non-commutative gauge theory in coordinate space (3.24) is not immediately suit- 
able for computer simulations. Therefore we make use of the finite dimensional 
representation introduced at the end of the last Section. The model we will arrive at 
is the twisted Eguchi-Kawai model (TEK). We start with a discussion of the origins 
of this model. 

4.1 The twisted Eguchi-Kawai model 
4.1.1 History of the TEK 

In 1982 Eguchi and Kawai conjectured that standard (commutative) U{N) and 
SU (N) lattice gauge theory in the large N limit is equivalent to their dimensional 
reduction to d — (one point) [77]. The link variables are replaced by Ufj_{x) — > Ufj_, 
and the partition function of the Eguchi-Kawai model (EK) simplifies to 

f fldU^e-'-^^''^ 
^^=1 (4.1) 
Sek[U] = -/3 5^ Tr {U^ Ul Ul) , 

where /.i,u = 1 . . . d. Eguchi and Kawai proved that in the large limit both models 
have the same Schwingcr-Dyson equations — and therefore the same Wilson loops 
— if the global U{lY symmetry of the phases is not spontaneously broken, lnd — 2 
this symmetry is unbroken and the equivalence holds exactly, but in > 2 this is 
not the case anymore. There the equivalence holds only at strong coupling [78]. 

A way out of this restriction to two dimensions is using twisted boundary condi- 
tions rather than periodic ones [79] . It is known that the behavior of the partition 
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function at weak coupling differs significantly between twisted and periodic bounda- 
ries. Therefore one can expect that the ^7(1)'^ symmetry breaking at weak coupling 
could be cured. 

To obtain the twisted Eguchi-Kawai model one applies the transformation 

U^.{x) ^ f/;(x) = Z^U^{x) , e Z^^ (4.2) 

to the partition function with the standard Wilson gauge action [72] 

j PC/exp(-5[C/]) 

S[U] = -/3 ^ Tr {U^{x)U,{x + aji)Ul{x + au)U,{x)^) . 



(4.3) 



The can in general depend on space-time, but for the purpose here it is enough 
to keep just the dependence on the orientation. The integration measure is invariant 
under this transformation and the transformed action reads 

SV\ --PYI ^/^-Tr {U^{x)U,{x + ail)Ulix + ai>)U,{x)^) . (4.4) 

The twist Zf^^, is the product of the factors Z^ at the boundary of the plaquette. It 
can be characterized by an integer valued anti-symmetric d x d matrix n^,^ 

2^^^^2.in,./N ^2:^ With n^. eZ. (4.5) 

Neglecting the space-time dependence of the hnk variables 1/^ then defines the TEK 
action 

and a general Wilson loop spanned by / x J plaquettes is given by 

W,.{I X J) = Z'jTr {Ul Ui Ul' Uy) . (4.7) 

In Ref. [80] it is shown that with this action the equivalence between the TEK and 
commutative U (N) and SU (N) lattice gauge theory also holds in > 2. 

Here wc are concerned with the 2d TEK. It seemed for a long time that adding 
a twist is not highly motivated in d = 2, since there even the ordinary EK model 
coincides with lattice gauge theory in the planar large N limit, which was solved 
analytically by Gross and Witten [81]. 

However, the situation changed suddenly due to a new interpretation of the TEK 
as an equivalent description of non-commutative gauge theory [82]. This equiva- 
lence was established by embedding the (dynamically generated) coordinates and 
momenta of the reduced model into matrices. These matrices can be mapped on 
functions, where the trace turns into an integral and the star products arise, so that 
one arrives at action (2.49). The construction of non-commutative U{n) gauge the- 
ories (for certain nG{l,2,...}) works out in this way at iV = cxd. ^ At finite N the 
conditions cannot be matched at the boundaries. 



^Note that N is here the size of the matrices 11^, whereas n refers to the gauge group U{n). 
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4.1.2 TEK at finite N 



An interpretation of the TEK at finite N occurred when Refs. [25, 73, 74] pointed 
out that the twisted reduced model can be mapped exactly on a non-commutative 
U (n) lattice gauge theory ^ in general dimension and rank of the gauge group. Here 
we will describe the case oi d — 2 and the star-unitary gauge group U{1). 

To show the equivalence we will use the finite dimensional representation of a 
non-commutative geometry, developed in Section 3.3. The starting point is the 
action of 2d non-commutative lattice gauge theory in terms of operators (3.24) 



Ui[D,U2D[ 



(4.8) 



Inserting the finite dimensional representation (3.34) of the general shift operator 



leads to 



,5 = -/3Tr 



The twist eaters obey the Weyl-'t Hooft commutation relation 

rir2 = Z*2V2T\ . 



(4.9) 



(4.10) 



The substitution (J^T ^ = in equation (4.9) along with the commutation relation 
(4.10) for the inner T's leads to the TEK action (4.6). 

In this representation the twist reads Z12 — — exp {7ri{N + l)/N). Compar- 
ing this with the twist needed for the TEK, where the exponent is given by 27ri/N 
multiplied by some integer, we find that (A^^-l- 1)/2 £ Z has to be fulfilled. Therefore 
we have to use odd values of N. 



Note that this interpretation of the TEK differs significantly from the original 
one. While in the original interpretation the space-time degrees of freedom become 

irrelevant in the large N limit, in the interpretation here the space-time degrees of 
freedom are exactly mapped onto matrix degrees of freedom already at finite N. 



4.1.3 Continuum limits 

The continuum limit a — also requires the limit N ^ 00, where the order of these 
limits is important for the resulting continuum theory. 

The plaricir limit 

In the planar limit one first sends N ^ 00, keeping the lattice spacing a fixed, 
followed by the continuum limit a — > along with /5 — cxd in a particular way 
dictated by the coupling constant renormalization. Equation (3.30) implies that in 
this limit ^ ^ 00 and only planar diagrams remain. In d = 2 there exists an exact 



^Remember that SU(n) gauge theories so far cannot be constructed on a non-commutative 
plane and therefore this equivalence exists only for U{n) gauge theories. 
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solution of planar U{oo) and SU{oo) [81]: rectangular Wilson loops follow the area 
law 

w{I X J) = exp(-K(/3)/J) , (4.11) 

where w{I x J) is the expectation value of the Wilson loops. The string tension 
k{(3) is given by 

<P) = < 1^ .-.I ■ (4-12) 



-ln(l-^) /3>i 



Equation (4.11) shows how the bare coupling has to be tuned as a function of a 
when one takes the continuum limit. In this case the scaling is exact if one identifies 

the lattice spacing a as 

a^^M0). (4.13) 

The equivalence between commutative U{N) and SU{N) theory and the TEK 
(and in two dimensions also EK) states that 

lim l^iWuil X J)) = w{I X J) . (4.14) 

This agreement was already shown in Ref. [77,83] and later on in Ref. [84] were finite 
effects where studied. 



The double scaling limit 

In the double scaling limit one takes the thermodynamic limit iV — oo and the 
continuum limit a — > at the same time keeping Na'^ = ttO fixed. Here we 
used the relation (3.30) between 9 and the lattice spacing. This leads to a fi- 
nite non-commutativity parameter 9 in the continuum limit and one recovers non- 
commutative 2d U{1) theory in M.'^. 

This is the continuum limit we studied. f3 has to be scaled as a function of the 
lattice spacing a in exactly the same way as in the planar theory. From equation 
(4.12) we read off /3 oc at large /3. Hence in d = 2 we are going to search for 
a double scaling limit keeping the ratio N/P constant. The question of renormaliz- 
ability of 2d non-commutative gauge theory at finite 9 can be answered by studying 
this large N limit of the TEK. If the observables converge to finite values, we can 
conclude that 2d lattice non-commutative gauge theory does have a finite continuum 
limit. Our results will be presented in the next Section. 



4.2 2d non— commutative U{1) theory 
4.2.1 The model 

We study the double scaling limit in two dimensional non-commutative C/(l) lattice 
theory. The explicit action in terms of dimensionless quantities reads 

Stek[U] = -NP J2 ^M-Tr {U^ Ul Ul) , (4.15) 
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with the twist given by (3.33). If not stated differently we set N//3 — 32. For details 
of the simulation see Appendix A. 2. 

Our main interest in this model is whether a finite continuum limit exists or not. 
Thus we study the double scaling limit — as described in the previous Section — of 
various observables. If we find a large iV scaling of these observables (at fixed Nj j3) 
this corresponds to a finite continuum limit. 

The observables we measure are constructed from Wilson loops spanned by / x J 
plaquettes 

W,,{I X J) = Z'^^Tr {Ul [// Ul' Uy) (4.16) 
and from open Polyakov lines, 

P,(I) = lY (U^) and P_^(/) = IV (u!!^) , (4.17) 

where / is the number of multiplied link variables and —/j, indicates the backward 
direction. 

In Section 3.2 we showed that one can construct star-gauge invariant observables 
from general open lines. Here we restrict ourselves to straight lines which are not 
necessarily closed by the boundary. 

4.2.2 Wilson loops and area law 

First we study the expectation values of Wilson loops. The expectation value of 
open Polyakov lines is vanishing. In fact the expectation value of any open contour 
vanishes due to the U{lY symmetry of the TEK. But n-point functions of these 
contours are sensible observables. 

In the EK model it was observed that square shaped Wilson loops converge faster 
to the known exact result than other rectangles with the same area [84]. Hence 
we also focus on square shaped Wilson loops, which corresponds to / = J in the 
definition (4.16). We define the normalized Wilson loop 

W{I) = j^{W,2{IxI))eC. (4.18) 

Due to the presence of the twist, there is no invariance under the transformation 

Ul ^ U2, U2 —>■ Ul. As a consequence Wf^u is in general complex and W12 = W21, 
hence the Wilson loop depends on the orientation. The real part represents the 
average over both orientations. Figure 4.1 shows the real part of the Wilson loops 
Re[iy(/)] as a function of the physical loop area. 

Large N scaling is clearly confirmed. At small up to moderate areas the Wilson 
loop follows the area law of planar lattice gauge theory given in Eq. (4.11), but at 
large areas it deviates and the real part oscillates around zero instead. This observa- 
tion is consistent over a wide range of N. Remarkably, not even the absolute value 
decays monotonously; at large areas it seems to fluctuate around an approximately 
constant value. This is illustrated in Figure 4.2 (top). 
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Figure 4.1: The real part of square shaped Wilson loops defined in Eq. (4.16) against 

the physical area. The plot at the top shows the Wilson loops for small areas, where 
the data follow the Gross- Witten area law (solid line). For large areas the real part 
of the Wilson loops deviate from this law and oscillates around zero instead. 
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Figure 4.2: Wilson loops against the physical area in polar coordinates. The solid 
line in the plot at the top shows the Gross-Witten area law. The phase of the Wilson 
loops increases linearly in the area beyond the area law regime (bottom). 
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Figure 4.2 (below) shows that beyond the Gross-Witten regime, the phase in- 
creases hnearly in the area. Additional measurements at N/j3 = 16,24,48, corre- 
sponding to different values of 9, are shown in Figure 4.3. They reveal that the 
phase $ of the Wilson loop W{I) is, in fact, given the simple relation 

$ = -^. (4.19) 

where the area is given by (a/)^. This relation holds to a high accuracy at large 
areas (a/)^ > 0{d), i.e. beyond the Gross-Witten regime. Relation (4.19) has been 
confirmed also for other (non-square) rectangular Wilson loops, which shows that 
the effect does not depend on the shape of the Wilson loops. ^ Indeed the formula 
(4.19) agrees with a Aharonov-Bohm effect in the presence of a constant magnetic 
field B — 1/9 across the plane. This is reminiscent of the description of non- 
commutative gauge theory by Seiberg and Witten [8]. This interpretation is also 
known in solid state physics. Our numerical results seem to support a picture of this 
kind. 

The observed large area behavior of the Wilson loops confirms that the continuum 
limit of non-commutative gauge theory is different from any ordinary gauge theory, 
hence we have found a new universality class. Since the non-locality in this model 
is of order v^, one might naively think that in the large area regime, area ^ \f9, 
the effect of non-commutativity is invisible. The fact that we the observe contrary 
can be understood as a manifestation of non-perturbative UV/IR mixing. 

In the two dimensional model there is no perturbative UV/IR mixing, since there 
are no UV divergent diagrams in the commutative case. Here we found this mixing 
at a completely non-perturbative level. 



In the double scaling limit of the untwisted Eguchi-Kawai model, the expectation 
values of the Wilson loops is real, and it remains positive even at large areas [84]. 
This means that the two models, twisted and untwisted, yield qualitatively different 
double scaling limits, although they become identical in the planar large N limit. 



4.2.3 2— point functions 

Here we consider 2-point functions. Figure 4.4 (top) shows the connected Wilson 
loop 2-point function 

(/) = {Wi2{I X I)W2i{I X /)) - {Wi2{I X I)){W2i{I X /)) e R, (4.20) 

again plotted against the physical area, for / = 1 . . . A^. In contrast to the Wilson 
loop itself, G^(/) is a real quantity, since both orientations of the Wilson loop are 
involved. 

In Figure 4.4 (below) we include a wave function renormalization factor 

q{W) ^ ^-0.6^(M/) _ ^4 21) 

^Note, however, that the expectation values of Wilson loops with the same area but with 
different shapes have in general different absolute values in the double scaling limit. 
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Figure 4.3: The phase of Wilson loops against the physical area for different values 
of 9. The phase increases linear in the area beyond the Gross-Witten regime for all 

e. 

The exponent —0.6 was found to be optimal for G^^^ to scale. Indeed it leads to a 
neat large scaling over more than two orders of magnitude in the physical area. 

Next we consider the Polyakov line 

P,{I) = Tr(u^), (4.22) 

which is also U{N) invariant and therefore has an interpretation as a star gauge 
invariant observable in non-commutative gauge theory. Their momentum p is related 
to the separation vector v between the two ends of the line. In general the relation 
is given by = Qf^^Pu modulo the periodicity of the torus, as discussed in Section 
(3.2). In the present case, the Polyakov line Pij,{I) corresponds to a momentum 
mode with p^, — 27ri/{Na), where the integer i is given by 7/2 and by {I + N)/2 
for even and odd /, respectively. In the following, we plot our results against the 
physical distance al for even I = 2, . . . , N — 1. 

The phase symmetry ^ makes (P^(/)) vanish, but the connected n-point functions 
(n > 1) of Polyakov lines are sensible observables. In Figure 4.5 we show the 2-point 
function ^ 

^In the terminology of non-commutative gauge theory, this corresponds to momentum conser- 
vation. 

^The choice of the direction fj, is irrelevant. In practice we average over both possibilities in 
order to enhance the statistics. 
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Figure 4.4: Two-point functions of square shaped Wilson loops (4.20) against the 
physical area. The raw data at the top do not show any scaling behavior, but with 
a wave function renormalization — > f3~^'^GY there is a clear scaling regime 
(bottom). 
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G^{I)^{P,{I)P.,{I))eR. 



(4.23) 



Note that there is no disconnected part in • Again we insert the wave function 
renormahzation which was optimal for the Wilson 2-point function 



As a function of the physical length al, the result is consistent with large scal- 
ing, as well as a universal wave function renormahzation. A similar wave function 
renormahzation was also observed in the EK model [84], where the optimal factor 
in relation (4.24) is modified to 

The large scaling of the Wilson loops as well as the large scaling of 2-point 
functions of Wilson loops and Polyakov lines, described in this Chapter, correspond 
to a finite continuum limit in 2d non-commutative U{1) gauge theory. This ob- 
servation therefore demonstrates the renormalizability of 2d non-commutative U{1) 
theory. 




(4.24) 
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Figure 4.5: Two-point functions of open Polyakov (4.23) against the physical area. 
The plot at the top shows that the raw data do not scale in N. Applying a wave 
function renormalization — > /3~°'^Gf leads to large A'' scaling in a significant 
interval. 



Chapter 5 

Numerical studies of the model 



The second model we investigated is the 3d non-commutative Xcj)^ theory. In Section 
2.2 we described the effects of the UV/IR mixing effect based on results of one loop 
calculations. We studied this model non-perturbatively and in this Chapter we 
present our results. For details of the simulations we refer to Appendix A. 3. 



5.1 Dimensionally reduced model 

Since we are in odd dimensions we cannot apply directly the construction of non- 
commutative field theories as described in Chapter 2 and 3. There the anti-symmetric 
non-commutativity tensor O had to be invertible. This restricts the dimension to 

be even. 

In addition to this rather technical problem related to odd dimensions, a non- 
commutative time could give rise to unitarity problems [85,86]. Ref. [87] showed that 
field theories on a non-commutative space satisfy the generalized unitarity relations. 
If the time is also non-commutative this is not the case. Therefore we exclude the 
time direction from non-commutativity. ^ 

Then the non-commutativity tensor is two dimensional and acts only in the 
2d non-commutative subspace. The star-product defined in (2.16) then reads 

dcf / ^ *' — — ^ \ 

(pi{x) -k 02(a;) = (pi{x) exp ( -diQijdj j (p2{x) with i,j = 1,2, (5.1) 

and the lattice version of the star-product is analogous. Here (f){x) means (f){xi,X2, t), 
where xi,X2 satisfy the commutation relation (2.18) and t commutes with all coor- 
dinates. 

The corresponding Weyl operators then also depend on the time (f){t) and the 
lattice action of this version of non-commutative theory reads, in analogy to 

^In four dimensions the problem is often avoided by taking two commuting and two non- 
commuting directions. 
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equation (3.19), 



t=l 



m 



A 



(5.2) 



where i = 1,2. There are now two kinetic terms: the first one uses the shift operators 
(3.11) to perform spatial translations in units of the lattice spacing a. The second 
kinetic term is the square of the standard discrete derivative in time direction. 

As in the TEK model we use here the finite dimensional representation (3.32). In 
this representation the Hermitian operators (f){t) turn into Hermitian matrices and 
the shift operators are replaced by the twist eaters Fj. 

Effectively we are mapping here a non-commutative A^"^ lattice theory, defined 
on a three dimensional N'^T lattice, to a one dimensional lattice with T sites. 







t 



t+1 



T 



Figure 5.1: The dimensionally reduced lattice. The fields (t){x\,X2,t) are mapped 
into Hermitian N x N matrices (j>{t). 



On each site there is a Hermitian N x N matrix (j){t) representing the Hermitian 
operators, see Figure 5.1. We use periodic boundary conditions 0(T) = 0(0). In our 
simulations we always set T — N. 

The dimensionless parameters and A in the action (5.2) can be identified with 
physical parameters 



2 2 2 ^ ^ 2 

m = a mphy, = — mp^y. 



A — aAphys — l/ ^ -^phys 



(5.3) 



Here we used the relation between the non-commutativity parameter 6 and the 
lattice spacing 9 = -a?N (3.30). With this identification the double scaling limit 
described in Subsection 4.1.3 leads to the non-commutative A0*^ model on M?. In this 
procedure the limits — > oo and a — > are taken such that a^N is kept constant, 
leading to a finite non-commutativity parameter. 

Here we will not study the continuum limit and the question of renormalizability. 
Instead we study the phase diagram and the UV/IR mixing effects in the regularized 
theory. 
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In Section 2.2.2 we discussed this issue in d = 4. Since the results depend on 
the dimension, we repeat the considerations of Section 2.2.2 for the case of three 
dimensions. 

We start with the non-planar one loop contribution to the irreducible two-point 
function in equation (2.44). In d — 3 this term reads 

rl» = + "V^) . (5.4) 

where A is again a momentum cut-off and p = {po,p)- Since we chose the time 
as commutative, only the spatial components p of the momentum p appear here. 
Introducing the effective cut-off AeS as in equation (2.38) 

and evaluating the Bessel function K1/2 leads the one loop corrected two-point 
function 

r(p) = p2 + + e A Aeff e"^ + O(A^) , (5.6) 

with ^ = 6(27r)3/^ • ^® absorbed the planar contribution into the effective 

mass Mcff. After removing the cut-off A in equation (5.6) we obtain the leading IR 
divergence 

T{p)^p' + M^^ + ^^^. (5.7) 

Again the two-point function is singular at zero momentum. However, here r(p) is 
not a function of p^, which leads to a different IR behavior of the theory. 



5.2 The phase diagram 

As we discussed in Section 2.2.3 the phase diagram of the non-commutative 3d \(f)'^ 
theory is expected to differ significantly from the phase diagram in the commutative 
case. In this Section we present our Monte Carlo results for this phase diagram. 



5.2.1 The order parameter 

When studying a phase diagram one first has to identify a suitable order parameter 
that indicates the symmetry breaking. In the model here we expect an Ising type 
phase, where the discrete symmetry 0(x) — > —4>{x) is broken spontaneously. In 
addition we expect — for sufficiently large coupling A — a striped phase, where 
the translation symmetry is broken spontaneously (see Section 2.2.3). Therefore we 
need an order parameter that is sensitive to both variants of symmetry breaking, to 
distinguish the two types of ordered phases. 
The momentum dependent quantity 

(5.8) 



M(m) = 



NT 



^0(m,i) 
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turned out to be a good choice. The vector m is the integer representation of the 
momenta introduced in equation (3.6). Here 0(m, t) is the spatial Fourier transform 
of the field 0(a;, t), where only the non-commutating coordinates are transformed. 
The expectation value of M{m) is zero in the disordered phase, where both symme- 
tries under consideration are unbroken. In an Ising type or uniform phase only the 
expectation value of M(0) is non-zero, since 



M(0)) 




^0(m = O,t) 



1 



\ iV2T 



t.x 



(5.9) 



is the standard order parameter of the spontaneous breakdown of the Z2 symmetry. 

— * 

A non-vanishing order parameter at m 7^ indicates a spontaneous breakdown 
of the translation symmetry and therefore it implies the striped or non-uniformly 
ordered phase. 

Since we are simulating the dimensionally reduced model in terms of the matrices 
0(i) we have to express this order parameter by these matrices. To this end we use 
the map defined in equation (3.38) and we obtain 



M(m) = 



NT 



^0(m,t) 



NT 



(5.10) 



This order parameter detects the disordered as well as the uniformly ordered 
phase. In the striped phase there will arise problems when the pattern for different 
configurations are rotated. For example if there are stripes for one configuration 
parallel to the Xi and for another configuration parallel to the X2 there would be 
two different non-vanishing order parameters. This effect is related to the rotation 
symmetry, which is also broken spontaneously in the striped phase. ^ To avoid this 
problem we defined the rotation invariant order parameter as the expectation value 
of 

M(A;) = max M(m) . (5.11) 

\rh\=k 

This order parameter depends only on the absolute value of the momentum and 
therefore the above mentioned problem does not occur. Based on this order param- 
eter we explored the phase diagram of the 3d A^"^ model. 



5.2.2 Numerical results 

With the algorithm described in Appendix A. 3 we generated configurations at var- 
ious values oi N, X and and measured the order parameter (5.11). The result 
is the phase diagram plotted in Figure 5.2. The points connected by lines display 
the phase transition between disordered phase and the ordered regime. The ordered 

^Actually rotation symmetry is explicitly broken on the lattice. However, in the commuta- 
tive case this symmetry is reduced to an invariance under rotations oi it/2, which can be broken 
spontaneously in the non-commutative lattice theory. 
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100 200 300 400 500 600 700 800 

Figure 5.2: The phase diagram of the non-commutative 3d A(/)^ in the - A plane. 

regime clearly splits into a uniform and in a striped phase, where the transition 
region between these two phases is marked by two vertical lines for each value of 
N. For each N the left line represents the largest value of A at which we are still 
in the uniform phase, and the right lines show the smallest value of A, where we are 
already in the striped phase. Both transitions stabilize in N if we multiply the axes 
by N^. 

To illustrate the different phases we mapped the matrices (j){t) back to coordinate 
space t)) using the map (3.34). We chose configurations in the four areas that 
can be distinguished. Some example snapshots of t) are shown in Figure 5.3. 




(a) N^m^^ -22.5 (h) N^m^ = -225 (c) N'^m^ = -90 (d) TV^m^ = -400 



Figure 5.3: Snapshots of single configurations (j){x,t) at some time t, at A'' = 45, 
iV^A = 90 (a,b) and = 450 (c,d). 

The dotted areas indicate t) > and in the blank areas t) is negative. 
Figure 5.3(a) and 5.3(c) show 0(x, t) in the disordered phase. The positive and 
negative spread all over as one expects in the disordered phase. Also in 

the uniformly ordered phase we find the expected behavior; 4>{x,t) is either positive 
for all X or negative for all x (Figure 5.3(b)). 
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Figure 5.3(d) shows a typical pattern with two stripes. In the range of the 
parameters plotted in the phase diagram 5.2 we always found two stripes (in the 
non-uniformly ordered phase). They were either parallel to the Xi or to the X2 axis. 
We will discuss the number of stripes separately in the next Subsection. In the rest of 
this Subsection we discuss the measurements that the phase diagram 5.2 is based on. 



To localize the phase transitions we started simulations in the disordered phase. 
After we reached equilibrium we measured M{k). Then we decreased slowly 
towards the expected phase transition. Figure 5.4 shows two example measurements 




Figure 5.4: The momentum dependent order parameter (M(fc)) against N'^m? at 
= 35. On the left we fixed at A^^A = 70, which leads to the uniform phase and on 
the right we are at N'^X = 350, leading to the striped phase. 



at = 35. We plotted in both cases the order parameter at = and k = 1 
against N'^m?. On the left we started in the disordered phase at N'^X = 70. Below a 
critical value of m^, we see a clearly non-vanishing {M{k)) at k = 0, while (M(l)) 
is zero for all values of m^. This indicates that we are entering the uniform phase. 
On the right we started in the disordered phase at A^^A = 350. Lowering leads 
to a non-zero order parameter at k = 1 and the standard order parameter (M(0)) 
is zero. The fact that (M(l)) ^ implies that the system is in the striped phase. 

Measurements of this kind allowed us to separate the uniform phase from the 
striped phase. ^ Whenever we see a dominant order parameter with k ^ we 
are in the striped phase. Unfortunately we were not able to determine an accurate 
transition line between these phases. We obtained a transition region, in which it 
was not possible to identify the order of the system. Depending on the starting 
conditions (see Appendix A) we found indication for both phases. We come back to 
this point in Section 5.4. 

From plots like shown in Figure 5.4 one can also identify the transition line 
between disordered phase and ordered regime, at least roughly. Whenever any order 
parameter becomes non-zero we hit the transition line. However, one can localize 
the transition more accurately by looking at the connected part of the two-point 
function of M{k), 

{M{kf)c = mkf) - (M(fc))^ . (5.12) 

•^A natural method the determine the transition line between uniform and striped phase is to 
start in the uniform phase and increase A or vice versa. However, it turned out that the initial 
patterns remain stable when the transition region is crossed. 
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From statistical mechanics it is known that this two-point function has a peak at the 
phase transition. Figure 5.5 shows two examples of these measurements, where we 
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Figure 5.5: The connected two-point function of M{k) defined in equation (5.12) 
against N'^m^ at N = 35. On the left at N'^X = 70 and on the right at N'^X = 350. 

used the same configurations as we did in Figure 5.4. This allowed us the determine 
the transition line between disordered phase and ordered regime to a high accuracy. 



5.2.3 The striped phase 

In Figure 5.4 we have seen that the two types of ordered phases can be distinguished 
by the momentum k of the non-vanishing order parameter {M{k)). We showed the 
two possibilities k = 0,1. It was left as an open question how {M{k)) behaves for 
other values of k. This we want to discuss here. 

In Figure 5.6 we plotted two examples of the order parameter against the mo- 
mentum k, where in both cases we are clearly in the ordered regime. On the left 
the system is in the uniformly ordered phase. Only the standard order parameter 
is non-zero, for all other momenta the order parameter is zero. This defines the 
uniform phase. 

The right plot in Figure 5.6 shows the system in the striped phase. There we 
see again a non-zero order parameter at = 1. In addition we also see a signal for 
(M(3)). This may appear as an indication for an underlying multi-stripe structure. 
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Figure 5.6: The momentum dependent order parameter (Af (fc)) against fc at = 
35. On the left in the uniform phase at iV^A = 70 and N^m^ = —210, and on the 
right in the striped phase at A^^A = 350 and iV^m^ = -400. 



which is not visible when we plot the sign of t), as we did in Figure 5.3. In fact 
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this is not the case. To illustrate this we define the function 

^i^i) = TLXl'^^^l'^2,^)• 



X2,t 



In this definition we assume that the stripes are parallel to the X2 axis. The function 
$(a;i) averages t) over X2 and t at Xi. If there is an additional structure ^ this 
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Figure 5.7: The momentum dependent order parameter {M{k)) against fc at = 35 
and <I>(xi) in the striped phase. On the left close to the phase transition at N'^X = 350 
and N'^m? = —250 and on the right at large negative values of m? in the striped 
phase at N'^X = 350 and N'^m^ = -620. 



should be seen in the behavior of To examine this we first computed the order 

parameter from configurations close to the phase transition and from configurations 
at large negative values of im?. The results are shown in Figure 5.7 (top). On the 
left the system is close to the phase transition and only (M(l)) is non-zero. On 
the right, at large negative values of rn?, we see that also (M(3)) and (M(5)) are 
different from zero. It is evident that when is further decreased, {M{k)) will also 
be non-vanishing for larger (odd) momenta k. 

We computed $(xi) on one configurations in each of these regions, to study how 
the additional momentum modes infiuence this function. The results are plotted 
in Figure 5.7 (bottom). If there is only one non-vanishing mode of the field this 
corresponds to a pure sine behavior. This is clearly confirmed by the plot on the 
left. The dashed line is a fit of the data to a sine function. ^ On the right ^{xi) 
is plotted in the case of various non-vanishing momentum modes. Here we see that 

''Of course the method will not find all possible patterns, but we also looked for additional 
structures in the other direction, which were not present. 

^The function $(a;i) is an average over X2 and t. The errors displayed in Figure 5.7 are 
computed with respect to this average with the binning method (see Appendix A.l). 
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these additional modes do not indicate a more complicated structure. It simply 
displays that deeply in the striped phase $(xi) is no more a pure sine, but the 
Fourier series involves higher momenta. As a fit function in this plot we used the 
first three odd terms in such an expansion, according to the three non-vanishing 
modes in Figure 5.7 (right, top), which fits the data perfectly. This result means 
for the structure of the field (^{x, t) that at large negative values of the number 
of stripes does not change, but the field is here restricted to only two values. At 
-C the kinetic term is negligible and the anti-ferromagnetic coupling dominates 
the system. 

As we mentioned above we obtained — in the range of the parameters shown 
in the phase diagram 5.2 — never more than two stripes. Analytically one expects 
not only two stripes, but also more complicated patterns [27]. To be in agreement 
with the perturbative conjecture these complex patterns should also occur on the 
lattice. The multi-stripe patterns are expected either on larger lattices or at larger 
non-commutativity parameter 6 or at larger coupling A. In fact for increased A 




(a) A = 10 (b) A = 10 (c) A = 100 (d) A = 100 



Figure 5.8: Snapshots of single configurations 4'{x, t) at a certain time t, at N — 35, 
A = 10 and = —6 (a,b) and at A = 100 and — —40. Tlie different patterns at 
eacli values of the parameters correspond to different starting configuration. 

we obtain also diagonal stripes, which support the picture of Gubser and Sondhi. 
However, only in the continuum limit, where the number of stripes should diverge, 
this picture can be ultimately verified. Some example snapshots at = 35 are 
displayed in Figure 5.8. The different patterns shown here correspond to different 
starting configurations. We postpone the discussion of the vacuum until we have 
considered the dispersion relation. 

5.3 Correlation functions 

In this Section we present our results concerning the correlation functions. The main 
interest here is to study the dispersion relation, but also the behavior of two-point 
functions in position space. 

5.3.1 Spatial correlators 

There are several predictions about the behavior of the correlation functions in 
coordinate space in the 4d A0'^ theory. It is expected from one loop perturbation 
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theory that the decay of this correlator at small A is not exponential [4]. For larger 
A the exponential behavior is restored but now depending on the coupling instead of 
the mass. This refers to the disordered phase. In the non-uniformly ordered phase 
an oscillation of the correlator is expected. Here we are studying the behavior of 
these correlators in d = 3. 

We are dealing here with one commuting and two non-commuting coordinates, 
namely the spatial coordinates. Since only the correlators in the non-commuting 
directions are expected to have an exotic behavior, we focus on spatial correlators. 
To be explicit, we study the correlation function ^ 



C{x) 



'^{(f){y,t)'k(j){y + x,t)) . 



(5.14) 



In Figure 5.9 the correlator (5.14) in the disordered phase at = 35 is shown. We 
measured the C{x) in one direction Xj keeping x fixed in the other direction. In 
Figure 5.9 we averaged over both directions. On the left we plotted the correlator 
(5.14) at A = 0.06 and on the right at A = 0.6. The decay is clearly not exponential 
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Figure 5.9: The correlator (5.14) against |x| in the disordered phase at A'' = 35. 
On the left the selfcoupling amounts to A = 0.06 and = —0.015. On the right the 
correlator is plotted at A = 0.6 and — —0.15. 



at these values of the coupling. 

In Figure 5.10 we show results at increased coupling. At moderately enlarged 
A the decay appears almost exponential (on the left), and at very large coupling 
the exponential decay is restored (on the right). The non-exponential decay of the 
spatial correlators was a result obtained from perturbation theory in four dimension. 
Here we observed this behavior in the 3d model. 



In the striped phase we have to take more care about the direction in which we 
measure the correlation function (5.14), since here we have patterns leading to a 
dependence of the correlator on the direction. Therefore an analysis of the stripe 
pattern has to be done first. This is achieved by evaluating the order parameter 

^In practice we use the twist eaters (3.32), which are a finite dimensional representation of the 
shift operators (3.11) to compute this quantity. We discuss this issue here and in the following in 
coordinate space since then the quantities are intuitively readable. 
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M{rh) defined in equation (5.8). This order parameter depends also on the orienta- 
tion of the condensed momentum mode and it is therefore not suitable for detecting 
the disordered - striped phase transition. However, here it is a perfect observable to 
detect the patterns. The momentum that maximizes M(m) dominates the pattern 
in the striped phase. For example, a maximum at m = (1, 0) indicates two stripes 
parallel to the X2 axis, and a maximum at m = (1, 1) or at m = (1, —1) indicates 
diagonal stripes. 

We measured the correlation function parallel and vertical to the stripes sepa- 
rately. For completeness we plotted in Figure 5.11 on the left also an example of 




1.2 

I 

0.8 
0.4 


-0.4 
-0.8 



parallel to the stripes r -G- -: 
vertical to the stripes ^ B--i 



a€»GoeeGeeoeeoQ€5eoeeeoeeo€X9eoe)ep 

\ 

B E3 

a 

EE3 







E3 







'B0bB' 



,0 



10 



15 _^ 20 

\x\ 



25 



30 



35 



Figure 5.11: The correlator (5.14) against \x\ in the ordered regime at = 35. On 
the left at A = 0.06 and = —0.1 in the uniform phase and on the right at A = 0.6 
and — —0.7 in the non-uniform phase. 



C(x) in the uniform phase. In this phase we see a strong correlation as expected in 
a uniform phase. However, in the striped phase (on the right) at A = 0.6 the spatial 
correlator behaves differently in the two directions. At this value of A we obtained 
two stripes parallel to one of the axis. Therefore we find in one direction (parallel 
to stripes) still a strong correlation. In the direction vertical to the stripes we see a 
strong anti-correlation. Since at this value of the coupling we have only two stripes, 
the correlator does not oscillate. 



The situation changes if we increase the coupling further. In Figure 5.12 C{x) 
is plotted at A = 10 and 100. In the last Section we showed already snapshots 
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of (j){x, t) at these values of the coupling. We saw there that depending on the 
starting configuration we obtain qualitatively different patterns. This leads to a 
different behavior of the spatial correlator (5.14). Again we refer to Section 5.4 for 
a discussion. Here we show results for the different striped phases separately. 

At these values of the coupling we find the maximum of M(m) at m = (1, 0) or 
at m = (1, 1) corresponding to two stripes parallel to the X2-axis or two diagonal 
stripes, respectively. In the first case the correlator C{x) behaves in the same way 
as in Figure 5.11 on the right. In the case of diagonal stripes we obtain a clear 
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Figure 5.12: The correlator (5.14) against \x\ in the striped phase at = 35. On 
the left at A = 10 and = -4 and on the right at A = 100 and = -40. The 
maximum of the order parameter is here at m = (1,1). 



oscillation of the correlation function vertical to the stripes as shown in Figure 5.12. 
Parallel to the stripes C{x) is strongly correlated. Note that here the oscillations do 
not indicate more than two stripes. 



5.3.2 Dispersion relation 

In Section 2.4 we discussed the broken Lorentz symmetry and the corresponding 
deformation of the dispersion relation in non-commutative theory, in the frame- 
work of perturbation theory. 

In a Lorentz invariant theory the energy squared is linear in p ^ and the dispersion 
relation is given by 

E{p)'=p' + M^s. (5.15) 
where Mgfj is the effective mass. In the non-commutative case equation (5.6) implies 
a dispersion relation of the form 

Eip)'=f + M^s + ^j^^, (5.16) 

where the last term is the leading one loop IR divergence. 

Here we study the dispersion relation non-perturbatively. Since we do not con- 
sider renormalization aspects in this model, we are only interested in the momentum 
dependence of the energy. To this end we considered the two-point correlation func- 
tion in time direction 

G'(m,r) = ^5^(Re(0*(m,t)0(m,t + r))) , (5.17) 
t 
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where 0(m, t) is again the spatial Fourier transform that we already used in the 
definition of the order parameter (5.8). On an infinite lattice G(m, r) decays expo- 
nentially in r. Since we are on a finite lattice with periodic boundary conditions 
(^(m, r + T) = (^(m, r), the correlator behaves like a cosh function, 



Gi^rh, t) oc (( 



-E{p)T 



+ e 



-Eip)iT-T) 



)• 



(5,18) 



Here E{p) is the energy, where we rescaled the integer representation of the momenta 
fh to their physical value p 

(5.19) 



P 



N 



-m . 



The correlator G(m, r) allows us to determine the momentum dependence of the 
energy, i.e. the dispersion relation, by studying its exponential decay. This can be 
done either by fitting G(m, r) to the function (5.18) or by studying the ratio of two 
subsequent values — log[G'(m, T + 1)/G{m, r)]. With the first method we obtain E{j>) 
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Figure 5.13: Determination of the energy E{p). On the left the data are fitted to 
the function (5.i8) and on the right we show — log(G(TO, t + 1)/G(m, t)). 



as a result of the fit, with the second method the energy is determined by a plateau. 
Results are shown in Figure 5.13. In this example the system is in the disordered 
phase and we computed G{m = 0, r). These measurements allow us to study the 
dispersion relation. 

We evaluated configurations close to the disordered - uniform transition and close 
to disordered - stripe transition. On all configurations we measured E{p) with the 
methods described above, for various momenta p. Two example results at = 45 
are displayed in Figure 5.14. On the left is linear in p^ as one expects in a Lorentz 
invariant theory. The solid line in this plot is the result of a fit to the dispersion 
relation (5.15), where we used the effective mass Mgfj as the only free parameter. 
Since we are on a finite lattice we see at larger momenta a deviation. Here the 
momentum dependence of the energy is E"^ oc (2sin|p|/2)^ (if Lorentz symmetry 
holds). A fit to this function is represented by the dashed line in Figure 5.14. This 
line fits the data very well. Since the energy minimum is at \p \ = at this value of 
A, we end up in the uniform phase when we decrease m^. 

The situation is different for larger A. This is shown in Figure 5.14 on the right. 
In the vicinity of |p | = we see here a clear deviation from the dispersion relation 
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Figure 5.14: The energy E'^{p) in the disordered phase at = 45, on the left at 
N'^X = 70 and N'^rn^ = -17.5, and on the right at N'^X = 700 and N'^m^ = -280. 



(5.15), at large momenta the linear behavior is restored. The increased energy at zero 
momentum is in full agreement with the perturbatively predicted IR divergence. Due 
to the discrete compactification (see Section 3.1) a finite lattice spacing serves also 
as an IR cut-off, along with the finite volume. Therefore we do not see a divergence. 
Since the cut-off influences the IR behavior we included an ad hoc cut-off as a 
consequence of equation (2.38). We fltted the data to the function 

E{pf = + + + a2{p + k) . (5.20) 

It turned out that the leading IR divergence in equation (5.16) does not describe 
the low momentum behavior sufliciently well. Therefore we added a linear term ^ 
in p from the expansion of the exponential function in equation (5.6). The relative 
errors of the flt parameters are shown in Table 5.1. The errors of the parameters 
of the deformation are rather large, which is not surprising since we have only one 
data point displaying this effect. 

However, the minimum of the energy is here clearly at the smallest non-vanishing 
momentum. This minimum implies that for decreased m? the system will be in the 
striped phase and the vacuum pattern will have two stripes parallel to one of the 
axes. For |p | = the energy increases again, which may indicate the IR divergence 
(at N oo). 

Here we observe that at small A, and equivalently for small 6, the effects of 
UV/IR mixing are strongly suppressed and we restore the behavior of the commuta- 
tive theory, including an Ising type phase transition. At larger A the UV/IR mixing 
effects become dominant and the phase transition changes its nature. This is in 
qualitatively agreement with the conjecture by Gubser and Sondhi [27]. 

For further increased A one expects the minimum at larger momenta. In addition 
E{0) should increase with A. Since at a slightly larger coupling the expected effect 
did not show up, we increased coupling drastically from A = 0.6 (the largest value 
at = 35 that is plotted in the phase diagram 5.2) to A = 10 and 100. The results 
are shown in Figure 5.15. We clearly see in these plots that E{0) increases with A. 

^Adding more terms in the expansion does not change the resufts within the errors. 
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Figure 5.15: The energy E'^{p) in tlie disordered pliase at = 35. On the left we 
are at A = 10 and = —1, and on the right at A = 100 and = —10. The energy 
at p = increases with A. 

We fitted the data again to the fit function (5.20). The results are the solid fines 
in Figures 5.15 and 5.16 and the relative errors of the fit parameters are shown in 
Table 5.1 in the last two lines. At these values of A also the errors of the parameters 
describing the deformation of the dispersion relation are under control. The results 
of the fits are consistent with the one loop result of perturbation theory (5.16) over a 
wide range of A. This is an unexpected result, since at these values of A effects from 
higher order perturbation theory are expected. One might conclude here that there 
are no qualitatively new IR singularities from higher loop contributions. However, 
this needs confirmation on larger lattices. 

To identify the minimum of the energy we plotted the dispersion relation in a 
smaller range (Figure 5.16). On the left, at A = 10, the values of the energy at 
k = \m\ = 1 and k = \m\ = -\/2 are equal within the statistical errors. From our 



1-5 I 1 1 1 1 1 1.7 




I I I I I I 2 I I I I I 

0.1 0.2 0.3 0.4 0.5 ' 0.1 0.2 0.3 0.4 0.5 



Figure 5.16: The energy E'^{p) in the disordered phase at iV = 35. On the left at 
A = 10 and = — 1, at fc = 1 and k = the energy is equal within the statistical 
errors. On the right at A = 100 and = —10, the minimum of the energy is shifted 
from A: = 1 to fc = V^. 

data it is not possible to predict which of these momentum modes will condense for 
decreased m^. On the right in Figure 5.16, at A = 100, the minimum of the energy 
is clearly at the second smallest non-vanishing momentum, which corresponds to 
k = \rh\ = \/2. Therefore the pattern in the non-uniform phase will have two 
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diagonal stripes. 



parameter 


2 

m 




ai 


a2 


relative error in Figure 5.14 right 


0.05 


0.25 


0.2 


0.2 


relative error in Figure 5.15 left 


0.02 


0.11 


0.08 


0.12 


relative error in Figure 5.15 right 


0.02 


0.15 


0.13 


0.15 



Table 5.1: The relative errors of the fit parameters. 



5.4 The phase diagram revisited 

As we have seen in Subsection 5.3.2, the investigation of the dispersion relation in 
the disordered phase provides information about the ordered regime. The mode 
that drives the phase transition, and therefore indicates the pattern of the ordered 
regime, is given by the momentum that minimizes the energy. 

This might give more insight into the transition region between uniform and 
striped phase. The difference between the energy at m = (0, 0) and at m = (1, 0) 
indicates the type of ordering that will occur in the ordered regime. In the uniform 
phase this difference is negative and in the striped phase positive. Equal energies 
indicate the phase transition. 

Since this is a small effect we performed high statistic ^ simulations at = 25 
in the range of A, which is marked as transition region in Figure 5.2. On these con- 
figurations we measured the energy gap E{0) —E{1). The result is plotted in Figure 
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Figure 5.17: The difference E{Q) — E{1) in the disordered phase at = 25 against 
N'^X. The vertical lines mark the transition region between uniform and striped 
phase, which was determined by measurements of the order parameter. 

5.17, where the region between vertical lines represent the aforementioned transition 
region. Outside this region and at its boundaries we clearly identify the uniform 
phase (left) and the striped phase (right). Inside the region the energy gap is zero 
within the errors with one exception close to the striped phase. With this method 
we cannot identify a transition line either. However, these measurements clarify 

^To extract this information we used 20000 configurations. 
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why here the pattern depends on the starting configurations. If E{0) and E{1) are 
almost equal it is (numerically) not clear which mode drives the phase transitions. 
Since numerical studies always suffer from the finite accuracy of computer numbers, 
it is hardly possible to resolve these small effects. Therefore already for two slightly 
different configurations the system might end up in one case in the uniform phase 
and in the other case in the striped phase, both of which appear to be stable. 

Regions where the minimum of the energy is difficult the identify are of course not 
restricted to the case discussed here. This occurs whenever the difference between 
the minimum of the energy and the next larger value of the energy is small compared 
to the absolute value of the energies. This is certainly the case when the minimum 
of the energy is changing with respect to the momentum. This obviously takes place 
in Figure 5.16 on the left. We also see this behavior at large A as Figure 5.16 on 
the right shows. In the first case we face the same problem as in the uniform - 
non-uniform transition area, which does not allow us to determine the behavior in 
the ordered regime. 

In the latter case the two lowest values of the energy are too close to resolve 
the correct pattern in the striped phase with our algorithm. However, from the 
analysis of the dispersion relation in the disordered phase we know that at = 35 
and A = 100 the striped phase has (1, l)-patterns. The additional (1, 0)-patterns in 
Figure 5.8 are only meta-stable. 

It is evident that more complex pattern will occur at further increased A. How- 
ever, in the limit of strong coupling the kinetic term in the action (5.2) becomes 
irrelevant. In Monte Carlo studies this causes a dramatic increase of the simulations 
steps that are needed to achieve the equilibrium (see Appendix A). 

Finally we comment on the orders of the phase transitions. We did not study this 
topic systematically, but we have some indications of which orders the transitions 
could be. 

The phase transition between disordered and ordered phase seems to be most 
likely of second order. We assume that because we do not see any indication of 
hysteresis at this transition. This holds for both, the disordered - uniform and 
the disordered - non-uniform transition. To study the hysteresis we performed 
simulations starting from the disordered phase and by slowly decreasing we 
entered the ordered regime. Once the system was clearly in the ordered regime 
we increased slowly towards the disordered phase. On these configurations we 
measured the order parameter. In a first order phase transition the order parameter 
would be different on both way. Since this hysteresis did not show up we assume a 
second order phase transition. 

We did not found a convincing indication of the order of the transition from the 
uniform to the striped phase. Here we obtained a transition region which prevents 
a prediction of the order. 
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Chapter 6 

The 2d non— commutative scalar 
model 



The occurrence of stripes in the ground state imphes the spontaneous breakdown 
of translation invariance. The case d = 2 is particularly interesting in this respect. 
Gubser and Sondhi argued based on an action of the Brazovskiian form [28] and the 
Mermin- Wagner theorem [29-31] that stripes cannot be stable in d = 2. However, 
Ambj0rn and Catterall pointed out that this theorem is not applicable, because 
here we deal with a non-local action as the star-product shows. In fact they did 
observe non-uniform patterns in their numerical results for non-commutative 
model ind — 2 [88], where the two coordinates obey the commutator relation (2.18). 
In general their results agree qualitatively with the results we got in d — 3. The 
only difference was that they obtained more complicated patterns than the two- 
stripe pattern at rather small coupling A. 

Inspired by their results we also studied the 2d version of the non-commutative 
model. In two dimensions the action reads 



where we used the twist-eaters (3.32) as shift operators We performed the 
similar measurements as we did in d — 3 beginning with the phase diagram. Again 
we used the order parameter (5.11), but since we are in two dimensions there is no 
sum over the time t. The phase diagram we obtained is plotted in Figure 6.1. As 
in the 3d case the ordered regime is split into a uniform and into a striped phase in 
agreement with Ambj0rn and Catterall. In this case it takes the factors N^^'^ and 
N'^ on the axes to stabihze the phase transitions in N. 

We also analyzed the striped phase in the range of the parameters that are plot- 
ted in Figure 6.1. In this range we always obtained patterns with two stripes. Since 
Ambj0rn and Catterall obtained multi-stripes already at these values of the cou- 
phng, we tried to understand the origins of these differences. According to Ref. [88] 




(6.1) 
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100 200 300 400 500 



Figure 6.1: The phase diagram of the non-commutative 2d Xcj)'^ in the - A plane. 

they used an algorithm that updates the complete matrix at once, in contrast to 
our algorithm where we updated pairs of matrix elements. Therefore we applied the 
same algorithm (although it has thermalization and ergodicity problems see Section 
A. 3) and indeed we also discovered other patterns than two stripes. Some exam- 
ple snapshots at = 35 are plotted in Figure 6.2. These patterns show up after 
approximately 500 simulation steps (starting form a random configuration) and sur- 
vived about 10^ further update steps. After increasing the number of update steps 
further to approximately 10^ to 10^ steps all patterns turn into the two stripe pat- 
tern (Figure 6.2 on the right) independent of the starting configuration. ^ We show 




Figure 6.2: Meta-stable patterns at iV = 35 and N'^X = 350 in the striped phase. 
After a very long thermahzation time all patterns turn into a two stripe pattern (on 
the right). 



here examples at A^ = 35, but the situation is the same at A^ = 45. However, the 
main result in Ref. [88] is of course the spontaneous breakdown of the translation 
invariance, and therefore the existence of a striped phase, which was not expected 

^In fact we faced the same problem in our study of the 3d model. In a first attempt we used 
the same algorithm as in Ref. [88] and we obtained multi stripes already at = 35. Since we 
saw a dependence of the patterns on the starting configuration we improved the algorithm, which 
unmasked these patterns as meta-stable. 
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in two dimensions. In this point our results fully agree with those of Ambj0rn and 
Catterall. 

We also computed the spatial correlation function defined in equation (5.14). As 
in d = 3 we obtained in the striped phase a strong correlation parallel to the stripes 
and a strong anti-correlation vertical to the stripes. 

For very large coupling A we finally obtained stable multi-stripe patterns. One 
example at = 45 is shown in Figure 6.3 on the left. For all starting configurations 




Figure 6.3: An example snapshot of a configuration on tlie left and the correlator 
(5.14) against \x\ on the right in the striped phase at = 45 and N'^X = 20000. 

we arrive at the same pattern, which persists for an apparently unlimited number 
of update steps even for the improved algorithm (A. 13). 

The pattern in Figure 6.3 has four stripes where the stripes themselves have 
some substructure. On the right we plotted the correlation function in position 
space (5.14). The stripes show up as an oscillation around zero of this correla- 
tor perpendicular to the stripes. The substructure of the stripes leads to a slight 
oscillation of the correlator parallel to the stripes. 
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Chapter 7 

Summary and conclusion 



We investigated non-commutative field theories in lower dimensions non-perturba- 
tively in the lattice approach. Since the non-local star-product arises also in the 
lattice formulation it is not suitable for a direct Monte Carlo investigation. There- 
fore we used a finite dimensional representation of the operator formulation of these 
field theories. This leads to dimensionally reduced models. 

In non-commutative gauge theory the reduced model (reduced to = 0) is 
given by the twisted Eguchi Kawai model (TEK). Within this model we investigated 
the continuum limit of 2d non-commutative U{1) theory. In the TEK this limit 
corresponds to the large N double scahng hmit. To this end we studied the double 
scaling limit of the 1-point function and the 2-point function of Wilson loops as well 
as the 2-point function of Polyakov lines. 

The first conclusion from our simulation results is that we do observe a dou- 
ble scaling limit as N, (3 oo. This corresponds to the continuum limit of the 
non-commutative U{1) gauge theory, which has therefore also finite observables. 
This observation demonstrates the non-perturbative renormalizability of 2d non- 
commutative U{1) gauge theory. 

The Wilson loop follows an area law at small physical areas, and in this regime 
non-commutative gauge theory agrees with planar standard gauge theory. However, 
at larger areas the Wilson loop becomes complex and the real part (the mean values 
over both loop orientations) begins to oscillate around zero. The phase is propor- 
tional to the physical area enclosed by the Wilson loop, irrespectively of its shape, 
and the coefficient of proportionality is given by the inverse of the non-commutativity 
parameter 6. This agrees with the Aharonov-Bohm effect in the presence of a con- 
stant magnetic background field 



Our results support this law, which is also a key element of the Seiberg and Witten 
description of non-commutative gauge theory [8]. Moreover, the same law also 
occurs in condensed matter physics: if one assumes a plane crossed by a constant 
magnetic flux, then the electrons in this plane can be projected to the lowest Landau 
level in a non-commutative space, where 9 — hc/eB [12-14]. 



65 



66 



Summciry and conclusion 



The behavior of the Wilson loop at large areas implies that we have found a 
qualitatively new universality class. At first sight, it may look surprising that the 
non-commutativity — which introduces a short-ranged non-locality in the action — 
changes the IR behavior of the gauge theory completely. However, this effect does 
not appear unnatural in the presence of UV/IR mixing. It is remarkable that this 
mixing effect does not occur in the perturbative expansion of this model, hence our 
results for large Wilson loops reveal a purely non-perturbative UV/IR mixing. 

For the connected Wilson loop 2-point function, as well as the 2-point function 
of the Polyakov hne, we can confirm the large N scahng. These observables are 
in agreement with a universal wave function renormalization, which yields a factor 
^-0.6 fQj. connected 2-point function. 

The second model we investigated was the 3d non-commutative theory, 
with a commutative time coordinate and two non-commutative space directions. 
Again the lattice model is mapped to a dimensionally reduced model, where here 
the reduced model is one dimensional. 

In this model we explored the phase diagram in the - A plane. In the ordered 
phase (at strongly negative m^) we found at small values of A a uniform order of 
the Ising type, as in the commutative case. At larger A — which amplifies the 
non-commutative effects — we observed striped patterns. Up to moderate values of 
A we obtained a pattern of two stripes parallel to the axes for = 15 ... 45. This 
corresponds to a minimum of the energy at the smallest non-zero lattice momentum. 
At very large A (compared to the coupling at the uniform - non-uniform phase 
transition) also other patterns showed up, in agreement with the conjecture of Gubser 
and Sondhi. However, the continuum limit has to be taken to confirm this agreement. 
We observed the same behavior in = 2 (see Chapter 6). Also there the multiple 
stripes occurred first at rather large values of A. The dominance of stripes implies 
the spontaneous breaking of translation symmetry, which is also possible in 2d since 
the action is non-local. 

In the ordered regime, the spatial correlations are dictated by the dominant pat- 
tern: uniform as in the commutative case, or striped with strong correlations in the 
direction of the stripes and strong anti-correlation vertical to them. This agrees with 
the prediction in Ref. [27]. In the disordered phase the spatial correlators deviate 
at small A from the exponential decay. At large A the decay is again exponential as 
was predicted from one loop perturbation theory in Ref. [4]. 

The correlators in momentum space do decay exponentially in time for all mo- 
menta. This property allowed us to study the dispersion relation in the disordered 
phase: at small A the dispersion relation behaves qualitatively like in the commu- 
tative case, but at large A there appears a jump at \p\ = as a non-commutative 
effect. Here we observed the energy minimum at \p \ — ^ for moderate A and at 
very large A the minimum is at \p \ = The results agree qualitatively with 

the results obtained in perturbation theory. Also from non-perturbative studies we 
observed an IR divergent behavior with the predicted l/\p\ divergence. However, 
since we obtained rather large fitting errors in the parameters of the IR dominant 
term this needs further confirmation. 
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Outlook 

In this work we presented results of our first steps in the non-perturbative study of 
non-commutative field theory. The next steps can be divided into next, near future 
and future projects. 

The next project is to study the 3d model at larger values of N. We expect 
from these studies more insight into the stripe structure in the ordered regime as 
well as into the IR behavior of the theory. Since the smallest non-zero momentum 
is given by |p| = a study at larger will resolve the interesting low momentum 
regime better. In addition we want to study the continuum limit and the related 
question of renormalizability. 

The near future project is to study non-commutative field theory in larger dimen- 
sions, for example 4d A(/)^ model or 3d gauge theory. Also a two or three dimensional 
a model would be interesting to study. In addition there are plans to study models 
including fermions for example a non-commutative Gross-Neveu model. 

In the long run we want to study models, which allow phenomenological predic- 
tions, like non-commutative pure gauge theory or non-commutative QED in d — A. 
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Appendix A 

The numerical methods 



In this appendix we discuss the details of the algorithms we used to simulate the 
non-commutative models. We start with some general comments on Monte Carlo 
simulations. 

A.l Monte Carlo simulations 

The idea of Monte Carlo (MC) simulations is to integrate approximately the infinite 
dimensional path integrals that appear in lattice formulation of any Euclidean field 
theory, with statistical methods. 

Vacuum expectation values of some observables 0[U] are computed in the path 
integral approach by 

(C») = ^ y DUO[U]e-^^^^ , (A.l) 

where S is the action. The action S may also depend on matter fields, ^ but for 
simplicity we restrict ourselves to gauge fields U^{x) E U{n). Z is the partition 
function 

Z^ J DUe-^^^^ . (A.2) 

A set of gauge fields {U}a, one for each link on the lattice, is a called a configu- 
ration. The idea of MC simulations is to generate as many configurations as possible 
via importance sampling. Importance sampling means that the configurations are 
generated according to the probability 

W[U] = ^e-^t^] . (A.3) 

Zj 

In addition to importance sampling a simulation algorithm has to fulfill the 
condition of ergodicity. This means that for any two configurations and {Uy^ 

the probability for a transition from one configuration out of the other, within finite 
number of update steps, is non-zero. 

^In that case also the matter fields have to be integrated. 
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On each configuration the observable 0[U] is calculated and the result is the 
value Oq:[?7]. The expectation value (A.l) is then given by 

1 ^ 

a=l 

In practice the number of configurations N is of course limited and therefore we get 
an approximation of the expectation value 

_ 1 ^ 

{0)^0 = -Y.O^[U], (A.5) 

a=l 

with an associated statistical uncertainty. If the configurations are statistically in- 
dependent this uncertainty reads 

i^Of = , (A.6) 

where O'^ is the average of O^. In practical applications the configurations gener- 
ated with a Monte Carlo simulation are always correlated and the error calculated 
with equation (A.6) is usually underestimated. This correlation in a sequence of 
configurations is called autocorrelation. For an infinitely large set of configurations 
the autocorrelation time is defined as 



'auto 



9 (O^) - (0)2 ■ V ■') 



There are several methods to take the autocorrelation time into account or to elimi- 
nate the effects of autocorrelation. Here we discuss only the methods that were used 
in this work. 



Error estimation 

In our evaluation of the sets of observables {Oa} we used two types of error estima- 
tions: the binning method and the jack-knife method. Both methods are based on 
dividing the set of observables into a number n of subsets. In each subset the average 
is calculated according to equation (A.5), where the number of configurations has 
to be replaced by the number of configurations in a subset. The result is a set of 
averages Oi, i = 1 . . .n. 

Within the binning method one considers the averages Oi as statistically inde- 
pendent and therefore the error is given by equation (A.6), where now the set of 
averages is Tiscd. The error calculated in this way varies with the number of bins n. 
In order not to underestimate the statistical error one should compute the errors for 
different values of n and take the largest error. 

The jack-knife method uses the averages Oi to construct n statistically indepen- 
dent subsets Ok, k — 1 . . .n in the following way. 
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Here Ok is the average of the Oi omitting the subset Ok- The error is now calculated 
according to equation (A. 6), using O^ as subset. This method has the advantage that 
it allows to compute the error of quantities that consist of two or more expectation 
values like the connected part of 2-point functions. 



Steps of a Monte Cetrlo simulation 

Before starting a simulation one has to perform tests of the algorithm. This is of 
course always important. Here it is indispensable, since in the case discussed there 
are no results available in the literature to compare qualitatively the obtained re- 
sults. The test for each algorithm will be discussed at the end of the corresponding 
Section. After these tests the simulations can be started. Here we briefly summarize 
the main steps of a typical Monte Carlo simulation. 



A Monte Carlo simulation consist of 

1. choosing a starting configuration. This can be either a random configuration 
or the gauge field set to certain values. 

2. performing update steps I: the first configurations are in general not distributed 
according to (A. 3), since a random configurations is located anywhere in con- 
figuration space. Therefore one performs update steps till the system is in 
equilibrium. These steps are called thermalization steps. 

3. performing update steps II: having reached the equilibrium a set of configura- 
tions is evaluated using equation (A. 5). The evaluated configurations should 
be separated by the autocorrelation time. 



A.2 Details of the TEK Model simulation 

For the TEK model we chose the heat hath algorithm [89]. In this method each link 
variable is updated separately, while the others are fixed. A new link is suggested 
and it is accepted with the probability 

P oc exp (-S[U]^ , (A.9) 

where S[U] is the part of the action that depends on U^. In other words, it is the 
sum over all plaquettes which involved U^, times the coupling /3. To enable this 
disentanglement the action has to be linear in the link variables. 

In the form of equation (4.6), the action cannot be simulated with the heat bath 
algorithm because it is not linear in the matrices U^. Following Ref. [90] we introduce 
an auxiliary field Q, which enters in a Gaussian form and linearizes the action in C/^, 



STEK[U,Q]^NpY,Tr 



H<i> 



(A.IO) 

{t;Mlui+t:^uiui)Q 
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The auxiliary field Q consists of general complex N x N matrices, and their integra- 
tion reproduces the TEK action (4.6), if t^^, = Z^^. The t/^ matrices are updated 
in SU{2) subgroups in the spirit of Cabbibo and Marinari [91], and the update of 
Qfj,i, is done by generating Gaussian variables. 

With this algorithm we generated configurations at various values of keeping 
N/f3 fixed. Table A.l shows the number of configurations we generated at each 
value of N. All configurations are separated by 50 update steps, which is more 



N 


25 


35 


55 


85 


125 


195 


255 


515 


# of configurations 


9000 


8400 


20000 


8000 


3400 


260 


120 


18 



Table A.l: Overview of the statistics in our simulation of the TEK model. 



than the autocorrelation time in these simulations (the maximum autocorrelation 
time measured was 40 update steps). We started the simulations with random 

configurations and it took approximately 100 update steps to reach equilibrium. 

The errors of the Wilson loops (4.18) and the Polyakov lines (4.22) were com- 
puted with the binning method. To estimate the errors of the two-point functions 
of Wilson loops (4.20) we used the jack-knife method, since this quantity consists 
of two expectation values. 

The simulations were performed on 10 Intel PCs with the operating system Linux. 
The PCs had a clock rate of 450 MHz. The generation of the configurations took 
approximately five month. 

For the measurement of the observables we had to compute high powers of 
matrices. Since they are not sparse matrices, in contrast to typical situations in 
lattice gauge theory, there is no simplified method to calculate these powers. One 
way to compute these powers is to evaluate the eigenvalues of the matrices first. All 
powers of the matrix are then just powers of the eigenvalues values. In this case one 
can use for example the conjugate gradient algorithm, which yields the eigenvalues 
up to a required precision. Unfortunately with this method we needed more matrix 
multiplications for a reasonable precision than with the direct multiplication. 

Therefore we decided the use a slightly modified direct multiplication of the ma- 
trices. The aim of the modification was to use more than one computer to calculate 
the powers of a given matrix without wasting too much computer time. ^ Assume 
that we want to compute f/^°°. To accelerate the calculation we want to use five 
computers, where to first computes to f/^° the second U^^ to f/^° and so on. 
This is of course only faster than using one computer if we manage to calculate the 
starting power for each computer in a fast way. 

For this we write the starting exponent in binary form. For example for the 
exponent 60 

60 = 1x2^ + 1x2^+1x2^ + 1x2^ + 0x2^ + 0x2°. (A.ll) 

In the next step these different powers of Uf^ are computed by taking repeatedly 
the square of the matrix. In the above example these are 5 matrix multiplications. 

^Note that we were not computing on a parallel machine. 
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The final result we get by multiplying the powers which have a non-vanishing co- 
efficient, here 4 multiplications, which means that we computed with 9 matrix 
multiplications. 

After the starting power of U/^ is generated the following powers are computed 
via direct multiplication. With this method it took approximately two months to 
evaluate the configurations in Table A.l. 

Tests of the algorithm 

In this model there exists a very stringent test of the implementation of the algo- 
rithm, namely the planar limit discussed in Section 4.1.3. In this limit there is an 
exact solution in the large N limit [81]. We reproduced these results in the planar 
limit with our code, confirming the validity of our investigations. In fact the pro- 
gram we used was only slightly modified compared to Ref. [84]. There the planar 
limit of the TEK model was studied, though, with a different twist and with even 
N, also reproducing the analytic solution of the planar limit. 

A. 3 Details of the A(/)^ model simulation 

For the scalar model we chose the standard Metropolis algorithm [92]. The idea of 
this algorithm is to suggest a new configuration, and accept this configuration with 
the probability 

P = min(l,exp(-(5„ew - S,,^))) . (A.12) 

If the action of the new configuration is smaller than the action of the old configura- 
tion it is accepted in any case, vice versa the new configuration is accepted with the 
probability exp(— AjS). The explicit update was performed by choosing an element 
of the Hermitian matrices (f){t) according to 

4>iM ^ 4>ij{t) = 4>ij{t) +er) and = . (A.13) 

Here 77 is a random complex number with (77) = and e e R"*" is an overall scale of 
the suggested changes, which has to be fixed. For every pair of complex conjugated 
elements we perform the accept-reject step (A.12). 

The parameter s should be chosen in a way to get a reasonable acceptance rate. If 
the acceptance rate is close to 1, the generated configurations are strongly correlated. 
If it is close to zero hardly any configurations are generated. A reasonable value 
would be 0.5, so one might try to fix £ to a value providing this acceptance rate. 

However, in this study we explore a phase diagram and the optimal value of 
e might depend on the different phases or it changes close to the phase transition. 
Therefore we decided to implement a dynamical value e: after updating the complete 
lattice we measure the acceptance rate. If the acceptance rate is above 0.6 we increase 
e by £/5, if the acceptance rate is below 0.3 we decrease e by e/5. ^ 



^Note that oscillations between two values of e are excluded with this choice. 
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from the interval 



y - ...... vc. L 0.5, 0.5] 

: A^~^/'*. In most cases e stabi- 



We choose the real and the imaginary part of 
times e. A reasonable starting value of e is then e 
lized after 3 — 10 update steps, depending on the coupling A, at a value close to the 
starting value. However, at very large values of the coupling e stabilized at a rather 
small value one order of magnitude smaller than N~^^^^ see below. 



The update of every independent matrix element of is of course more expen- 
sive than updating each matrix at once. However, it improves thermalization 
as well as ergodicity properties. In Figure A.l we show the history of the action for 
three different starting configurations. For the history shown by the thin line we 
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Figure A.l: The action is plotted against the number of update steps close to the 
disordered - ordered phase transition. The solid line shows the result when updating 
one complete matrix (j){t) at once, the thin line shows the result when updating every 
independent matrix element separately. 

used the update procedure (A. 13) and the thick lines represent the result when we 
update each matrix at once. In both methods we used a dynamical parameter e. 
First of all we see that the system thermalize much faster when updating the ele- 
ments separately. This already compensates the additional computer time needed 
for this algorithm. We also see that with the procedure, where we update one ma- 
trix at once, the action fluctuates much less than with the update procedure (A. 13). 
This means that we hardly proceed in configuration space and the configurations 
are strongly correlated. 

The reason for this is that if we change a complete matrix the change of the 
action is rather large and we get a reasonable acceptance rate only if we decrease 
the overall factor e to a tiny value. The configuration is hardly changed within an 
update step. 

In practice we also measured e as an technical parameter. It depends only slightly 
on w?, but it strongly depends on the coupling A. For very large couplings we have 
again an ergodicity problem. The values of A discussed in Chapter 5 can be consid- 
ered as the upper limit for the algorithm described here. 
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Besides the qualitative structure of the phase diagram, no details about the 
phases structure of this theory were known. As we expect a striped phase, with 
general patterns, we are possibly dealing in this phase with meta-stable states, 
which are usually difficult to identify. 

Therefore we always started from six independent (random) starting configu- 
rations and initialized the random number generator differently, i. e. we made for 
any set of A^, A, six independent simulations. This is of course much more ex- 
pensive than starting only from one configuration, because we have to thermalize 
all six runs. It has several benefits, however. If there are meta-stable states, they 
certainly depend on the starting conditions. With different starting conditions we 
should therefore see, in that case, different vacua. This could indicate that we should 
apply more thermalization steps. On the other hand, if we thermalize the system 
further, this shows that the algorithm is not ergodic, and we are stuck in a subset 
of configuration space. 

In fact, we saw these different meta-stable states, indicated by different patterns, 
at some stage of the simulations. At moderate values of A, applying more thermal- 
ization steps all starting conditions lead to the same results. At very large A this 
tunneling from meta-stable patterns to the stable pattern occurred only close to the 
ordered - disordered phase transition (see Section 5.4). 

A second advantage is that it clarifies the error estimation. As discussed in the 
last Section, the number of bins influences the statistical error in the binning as well 
as in the jack-knife method. In our case there is a natural choice of bins. Since we 
generate six statistically independent sets of configurations, we also chose six bins 
for the error analysis. 

Table A. 2 shows the number of configurations we generated with the algorithm 
described above. The necessary statistics to study the phase diagram is relatively 
small; this is shown in the second row of table A. 2. To compute the 2-point functions. 



N 




15 


25 


35 


45 


# of configurations 


phase diagram 


300 


180 


60 


30 


2-point functions 


10000 


10000 


5000 


3000 



Table A. 2: Overview of the statistics in our simulation of the 3d X<j)'^ model. 



especially those in momentum space, we needed significantly more configurations to 
obtain small statistical errors. At some values of the parameters we performed high 
statistic simulations where we doubled the number of configurations. The autocor- 
relation time in this simulation was in the range 10 - 90 update steps, depending on 
the observable and on the point in the phase diagram under consideration. Therefore 
we measured the observables from configurations separated by 100 update steps. 

It took approximately 10 month on 30 PCs (with a clock rate in the range 450 
- 1500 MHz) to generate the configurations in Table A. 2. The evaluation of the 
configuration did not increase the overall needed computer time, since we did not 
need high powers of the N xN matrices 4>{t). Most of the (usually expensive) matrix 
multiplications are in this evaluation multiplications with the constant twist eaters 
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(see Chapter 5). These multiphcation can be performed by shifting or modifying 
the phases of matrix elements of 4>, which is much faster than a direct multiphcation. 

Tests of the algorithm 

To test the algorithm is here more sophisticated than in the TEK model, since no 
exact solution is known in this case. Therefore we computed the expectation value 
of the action S in equation (5.2) to the first order in A , 

^ {S)^co{m) + Xci{m) + 0{X^). (A.14) 



For small A our algorithm should reproduce this result. We give now a brief account 
on how the coefficients cq and ci were obtained [93] . 

To compute cq we used the fact that the partition function is invariant under a 
change of variables 4>ij{t) — > (1 + €)4>ij{t) — (t>ij{t). The integration measure T>4> and 

the action S[(f)\ a=o transform under this shift like 

= (1 + e)^'^P0 and S[4>'] ^ {1 + efS[4>] , (A.15) 

where the transformation of the action is only valid in first order in A. Since the 
partition function is invariant under this transformation we obtain 



S[4>] _ / T)li ^-s[4,'] 



N^T J D0e-^t'^l - 2 J V(j)S[4)] e 



-s\4>\ 



+ O(e^) . 
(A.16) 



To satisfy this equation the term in square brackets has to vanish. This amounts to 
the condition 

^(5) = 0.5 = Co. (A.17) 
Note that the result for cq is independent of N ^ T and m^. 

For the determination of ci we extend the above consideration to the interacting 
case A 7^ 0. This leads to the identity 

Tyjorri (('Sfree) + 2(<S'interaction)) = 0.5 

^ ^ ^ (A.18) 

f^2rp + (-S'interaction)) = 0.5 . 

Therefore Ci is given by 



" ANT 



i^^5^Tr0(t)^^ . (A.19) 
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This expectation value can be calculated by transforming the partition function 
into momentum space, with the discrete Fourier transformation defined in equation 
(3.38) and using Wick's theorem to evaluate the resulting expectation value. The 
equation for the coefficient Ci is now given by a sum over momenta. The term that 
is summed over is very lengthy and complex, therefore we will not show the explicit 
result. However, for given parameters N, T and it can be computed for example 
with the program mathematica. 

To test if our algorithm provides this result, we computed the action at = 25 
and = 1. The coefficient cq is always given by cq = 0.5. For the chosen value 
of = 1 we find Ci = —0.0146. Figure A. 3 clearly confirms that our program 
reproduces this result. 




Figure A. 2: The expectation value of tlie action against A, at = 25 and — 1. 
The circles show the results of the simulation, the straight line is the expectation 
value action according to a first order expansion in A. 



78 Appendix A The numerical methods 



Bibliography 



[1] H. S. Snyder. Quantized Spacf^time. Phys. Rev. 71, 38 (1947). 

[2] C. N. Yang. On Quantized Space-time. Phys. Rev. 72, 874 (1947). 

[3] T. Filk. Divergencies in a Field Theory on Quantum Space. Phys. Lett. B376, 
53 (1996). 

[4] S. Minwalla, M. Van Raamsdonk, and N. Seiberg. Non-commutative Pertur- 
bative Dynamics. JHEP 02, 020 (2000), hep-th/9912072. 

[5] J. C. Varilly and J. M. Gracia-Bondia. On the Ultraviolet Behaviour of Quan- 
tum Fields over Non-conmiutative Manifolds. Int. J. Mod. Phys. A14, 1305 
(1999), hep-th/9804001. 

[6] A. Connes and M. Rieffel. Yang-Mills for Non-commutative Two-tori. Con- 
temp. Math. Oper. Algebra. Math. Phys. 62, 237 (1987). 

[7] A. Abouelsaood, C. Callan, C. R. Nappi, and S. A. Yost. Open Strings in 
Background Gauge Fields. Nucl. Phys. B280, 599 (1987). 

[8] N. Seiberg and E. Witten. String Theory and Non-commutative Geometry. 
JHEP 09, 032 (1999), hep-th/9908142. 

[9] S. Doplichcr, K. Fredenhagen, and J. E. Roberts. Space-time Quantization 
Induced by Classical Gravity. Phys. Lett. B331, 39 (1994). 

[10] S. Dophcher, K. Fredenhagen, and J. E. Roberts. The Quantum Structure of 
Space-time at the Planck Scale and Quantum Fields. Commun. Math. Phys. 
172, 187 (1995). 

[11] M. R. Douglas and N. A. Nekrasov. Non-commutative Field Theory. Rev. 
Mod. Phys. 73, 977 (2002), hep-th/0106048. 

[12] S. M. Girvin and A. H. MacDonald. Off Diagonal Long Range Order, Oblique 
Confinement, and the Fractional Quantum Hall Effect. Phys. Rev. Lett. 58, 
1252 (1987). 

[13] J. Bellissard, A. van Elst, and H. Schulz-Baldes. The Non-commutative Geom- 
etry of the Quantum Hall Effect. (1994), cond-mat/9411052. 



79 



80 



Bibliography 



[14] S. M. Girvin. The Quantum Hall Effect: novel Excitations and broken Symme- 
tries. (1998), cond-mat/9907002. 

[15] L. Susskind. The Quantum flail Fluid and Non-commutative Chern Simons 
Theory. (2001), hep-th/0101029. 

[16] A. P. Polychronakos. Quantum ffall States on the Cylinder as Unitary Matrix 
Chern-Simons Theory. JHEP 06, 070 (2001), hep-th/0106011. 

[17] W. Bietenholz, F. Hofheinz, and J. Nishimura. The Renormalizabihty of 2d 
Yang-Mills Theory on a Non-commutative Geometry. JHEP 09, 009 (2002), 
hep-th/0203151. 

[18] W. Bictculiolz. F. Hoflicinz. and J. Nishimura. Simulating Non-commutative 
Field Theory. Nucl. Phys. Proc. Suppl. 119, 941 (2003), hep-lat/0209021. 

[19] W. Bietenholz, F. Hofheinz, and J. Nishimura. Non-commutative Field The- 
ories beyond Perturbation Theory. Fortsch. Phys. 51, 745 (2003), hep- 
th/0212258. 

[20] K. Osterwalder and R. Schrader. Axioms for Euclidean Green's Functions. 2. 
Commun. Math. Phys. 42, 281 (1975). 

[21] D. Bahns, S. Dophcher, K. Fredenhagen, and G. PiaciteUi. On the Unitarity 
Problem in Space-time Non-commutative Theories. Phys. Lett. B533, 178 
(2002), hep-th/0201222. 

[22] C. Rim and J. Yee. Unitarity in Space-time Non-commutative Field Theories. 
(2002), hep-th/0205193. 

[23] Y. Liao and K. Sibold. Time-ordered Perturbation Theory on Non- 
commutative Space-time. II. Unitarity. Eur. Phys. J. C25, 479 (2002), hep- 
th/0206011. 

[24] A. Connes, M. R. Douglas, and A. Schwarz. Non-commutative Geometry 
and Matrix Theory: Compactification on Tori. JHEP 02, 003 (1998), hep- 
th/9711162. 

[25] J. Ambj0rn, Y. M. Makeenko, J. Nishimura, and R. J. Szabo. Lattice Gauge 
Fields and Discrete Non~commutative Yang-Mills Theory. JHEP 05, 023 
(2000), hep-th/0004147. 

[26] R. J. Szabo. Quantum Field Theory on Non-commutative Spaces. (2001), 
hep-th/0109162. 

[27] S. S. Gubser and S. L. Sondhi. Phase Structure of Non-commutative Scalar 
Field Theories. Nucl. Phys. B605, 395 (2001), hep-th/0006119. 

[28] S. Brazovkii. Phase Structure of an fsotropic System to a Non-uniform State. 
Zh. Eksp. Teor. Fiz A 15, 319 (1975). 



Bibliography 



81 



[29] N. D. Mermin and H. Wagner. Absence of Ferroinagnetism or Antiferromag- 
netism in One-dimensional or Two-dimensional Isotropic Heisenberg Models. 
Phys. Rev. Lett. 17, 1133 (1966). 

[30] P. C. Hohenberg. Existence of Long-Range Order in One and Two Dimensions. 
Phys. Rev. 158, 383 (1967). 

[31] S. R. Coleman. There are no Goldstone Bosons in Two Dimensions. Commun. 
Math. Phys. 31, 259 (1973). 

[32] G. Chen and Y. Wu. Renormalization Group Equations and the Lifshitz Point 
in Non-commutative Landau-Ginsburg Theory. Nucl. Phys. B622, 189 (2002), 
hep-th/0110134. 

[33] K. Matsubara. Restrictions on Gauge Groups in Non-commutative Gauge The- 
ory. Phys. Lett. B482, 417 (2000), hep-th/0003294. 

[34] A. Armoni. Comments on Perturbative Dynamics of Non-commutative Yang- 
Mills Theory. NucL Phys. B593, 229 (2001), hep-th/0005208. 

[35] M. Hayakawa. Perturbative Analysis on Infrared Aspects of Non-commutative 
QED on R**4. Phys. Lett. B478, 394 (2000), hep-th/9912094. 

[36] B. Jurco, S. Schraml, P. Schupp, and J. Wess. Enveloping Algebra Valued 
Gauge Transformations for Non-Abelian Gauge Groups on Non-commutative 
Spaces. Eur. Phys. J. C17, 521 (2000), hep-th/0006246. 

[37] C. Chu and H. Dorn. Non-commutative SU(N) and Gauge Invariant Baryon 
Operator. Phys. Lett. B524, 389 (2002), hep-th/0110147. 

[38] I. Chepelev and R. Roiban. Renormalization of Quantum Field Theories on 
Non-commutative R**d. T. Scalars. JHEP 05, 037 (2000), hep-th/9911098. 

[39] S. Sarkar. On the UV Renormahzability of Non-commutative Field theories. 
JHEP 06, 003 (2002), hep-th/0202171. 

[40] L. Griguolo and M. Pietroni. Wilsonian Renormahzation Group and the Non- 
commutative IR/UV Connection. JHEP 05, 032 (2001), hep-th/0104217. 

[41] L. Griguolo and M. Pietroni. Hard Non-commutative Loops Resummation. 
Phys. Rev. Lett. 88, 071601 (2002), hep-th/0102070. 

[42] H. O. Girotti, M. Gomes, V. O. Rivcllcs, and A. J. da Silva. A Consistent 
Non-commutative Field Theory: The Wess-Zumino Model. Nucl. Phys. B587, 
299 (2000), hep-th/0005272. 

[43] H. O. Girotti, M. Gomes, A. Yu. Petrov, V. O. Rivelles, and A. J. da Silva. The 
Three-Dimensional Non-commutative Nonlinear Sigma Model in Superspace. 
Phys. Lett. B521, 119 (2001), hep-th/0109222. 



82 



Bibliography 



[44] R. Wulkenhaar. Non-renorinalizability of 6-expanded Non-commutative QED. 
JHEP 03, 024 (2002), hep-th/0112248. 

[45] D. Colladay and V. A. Kostclccky. Lorcntz-violating Extension of the Standard 
Model. Phys. Rev. D58, 116002 (1998), hep-ph/9809521. 

[46] S. M. Carroll, J. A. Harvey, V. A. Kostclccky, C. D. Lane, and T. Okamoto. 
Non-commutative Field Theory and Lorentz Violation. Phys. Rev. Lett. 87, 
141601 (2001), hep-th/0105082. 

[47] G. Amclino-Camclia, L. Doplichcr, S. Nam, and Y. Sco. Phenomenology of 
particle production and propagation in string-motivated canonical noncommu- 
tative spacetime. (2001), hep-th/0109191. 

[48] A. Matusis, L. Susskind, and N. Toumbas. The IR/UV Connection in the 
Non-commutative Gauge Theories. JHEP 12, 002 (2000), hep-th/0002075. 

[49] G. Amelino-Camclia, J. Ellis, N. E. Mavromatos, D. V. Nanopoulos, and Subir 
Sarkar. Potential Sensitivity of Gamma-Ray Burster Observations to Wave 
Dispersion in Vacuo. Nature 393, 763 (1998), astro-ph/9712103. 

[50] S. D. Biller et al. Limits to Quantum Gravity Effects from Observations of TeV 
Flares in Active Galaxies. Phys. Rev. Lett. 83, 2108 (1999), gr-qc/9810044. 

[51] The HESS Project. http://www.mpi-hd.mpg.de/hfm/HESS/HESS.html. 

[52] T. Lohse, 02 2003. Private Communication. 

[53] A. de Angelis. GLAST, the gamma-ray large area space telescope. (2000), 
astro-ph/0009271. 

[54] F. Aharonian et al. TeV Gamma Rays from the Blazar H 1426+428 and the 
diffuse Extragalactic Background Radiation. (2002), astro-ph/0202072. 

[55] K. Greisen. End to the Cosmic Ray Spectrum? Phys. Rev. Lett. 16, 748 (1966). 

[56] F. W. Stecker and O. C. De Jager. New Upper Limits on Intergalactic Infrared 
Radiation from High-energy Astrophysics. Astrophys. J. 415, L71 (1993). 

[57] G. Zatsepin and V. Kuzmin. Upper Limit of the Spectrum of Cosmic Rays. 
JETP Lett. 4, 78 (1966). 

[58] N. Hayashida et al. Updated AGASA event list above 4*10**19-eV. Astrophys. 
J. 522, 225 (1999), astro-ph/0008102. 

[59] G. Amelino-Camelia and T. Piran. Planck-scale deformation of Lorentz sym- 
metry as a solution to the UHECR and the TeV-gamma paradoxes. Phys. Rev. 
D64, 036005 (2001), astro-ph/0008107. 

[60] J. L. Hewett, F. J. Petriello, and T. G. Rizzo. Signals for Non-commutative In- 
teractions at Linear Colliders. Phys. Rev. D64, 075012 (2001), hep-ph/0010354. 



Bibliography 



83 



[61] M. M. Sheikh- Jabbari. Renormalizability of the Supersyinmetric Yang-Mills 
Theories on the Non-commutative Torus. JHEP 06, 015 (1999), hep- 
th/9903107. 

[62] T. Krajewski and R. Wulkenhaar. Perturbativc Quantum Gauge Fields on 
the Non-commutative Torus. Int. J. Mod. Phys. A15, 1011 (2000), hep- 
th/9903187. 

[63] I. Ya. Aref'eva, D. M. Bclov, A. S. Koshelev, and O. A. Rytchkov. UV/IR 
Mixing for Non-commutative Complex Scalar Field Theory. 11: Interaction with 
Gauge Fields. Nucl. Phys. Proc. SuppL 102, 11 (2001), hep-th/0003176. 

[64] J. A. Aguilar-Saavedra et al. TESLA Technical Design Report Part 111: Physics 
at an e+e- Linear GoUider. (2001), hep-ph/0106315. 

[65] T. Abe et al. Linear GoUider Physics Resource Book for Snowmass 2001. 4: 
Theoretical, Accelerator, and Experimental Options. (2001), hep-ex/0 106058. 

[66] R. W. Assmann et aL A 3-TeVe+ e- Linear GoUider Based on GLIG Technology. 
SLAC-REPRINT-2000-096. 

[67] S. Godfrey and M. A. Doncheski. Signals for Non-commutative QED in 67 and 
-f-f GoUisions. Phys. Rev. D65, 015005 (2002), hep-ph/0108268. 

[68] T. G. Rizzo. Signals for Non-commutative QED at Future e+ e- Golliders. 
(2002), hep-ph/0203240. 

[69] I. Hinchhffe and N. Kersting. Review of the Phenomenology of Non- 
commutative Geometry. (2002), hep-ph/0205040. 

[70] A. Anisimov, T. Banks, M. Dine, and M. Graesser. Comments on 
Non-commutative Phenomenology. Phys. Rev. D65, 085032 (2002), hep- 
ph/0106356. 

[71] S. Godfrey and M. A. Doncheski. Phenomenology of Non-commutative Field 
Theories. (2002), hep-ph/02 11247. 

[72] K. G. Wilson. Confinement of Quarks. Phys. Rev. DIO, 2445 (1974). 

[73] J. Ambj0rn, Y. M. Makeenko, J. Nishimura, and R. J. Szabo. Finite N Ma- 
trix Models of Non-commutative Gauge Theory. JHEP 11, 029 (1999), hep- 
th/9911041. 

[74] J. Ambj0rn, Y. M. Makeenko, J. Nishimura, and R. J. Szabo. Non-perturbative 
Dynamics of Non-commutative Gauge Theory. Phys. Lett. B480, 399 (2000), 
hep-th/0002158. 

[75] R. J. Szabo. Discrete Non-commutative gauge theory. Mod. Phys. Lett. A16, 
367 (2001), hep-th/0101216. 



84 



Bibliography 



[76] P. van Baal and B. van Geemen. A Simple Construction of Twist Eating Solu- 
tions. J. Math. Phys. 27, 455 (1986). 

[77] T. Eguchi and H. Kawai. Reduction of Dynamical Degrees of Freedom in the 
Large N Gauge Theory. Phys. Rev. Lett. 48, 1063 (1982). 

[78] G. Bhanot, U. M. Hcllcr, and H. Neuberger. The Quenched Eguchi-Kawai 
Model. Phys. Lett. B113, 47 (1982). 

[79] A. Gonzalez-Arroyo and M. Okawa. A Twisted Model for Large N Lattice 
Gauge Theory. Phys. Lett. B120, 174 (1983). 

[80] A. Gonzalez-Arroyo and C. P. Korthals Altes. Reduced Model for Large N 
Continuum Field Theories. Phys. Lett. B131, 396 (1983). 

[81] D.J. Gross and E. Witten. Possible Third Order Phase Transition in the Large 
N Lattice Gauge Theory. Phys. Rev. D21, 446 (1980). 

[82] H. Aoki, N. Ishibashi, S. Iso, H. Kawai, Y Kitazawa, and T. Tada. Non- 
commutative Yang-Mills in IIB Matrix Model. Nucl. Phys. B565, 176 (2000), 
hep-th/9908141. 

[83] K. Fabricius and O. Haan. Numerical Analysis of the Twisted Eguchi-Kawai 
Model in Two-Dimensions. Phys. Lett. B131, 399 (1983). 

[84] T. Nakajima and J. Nishimura. Numerical Study of the Double Scaling Limit 
in Two- dimensional Large N Reduced Model. Nucl. Phys. B528, 355 (1998), 
hep-th/9802082. 

[85] N. Scibcrg, L. Susskind, and N. Toumbas. Space-time Non-commutativity and 
Causality. JHEP 06, 044 (2000), hep-th/0005015. 

[86] H. Bozkaya, P. Fischer, H. Grosse, M. Pitschmann, V. Putz, M. Schweda, and 
R. Wulkenhaar. Space-time Non-commutative Field Theories and Causality. 
(2002), hep-th/0209253. 

[87] J. Gomis and T. Mehen. Space-time Non-commutative Field Theories and 
Unitarity. Nucl. Phys. B591, 265 (2000), hep-th/0005129. 

[88] J. Ambj0rn and S. Cattcrall. Stripes from (Non-commutative) Stars. Phys. 
Lett. B549, 253 (2002), hep-lat/0209106. 

[89] M. Creutz. Monte Carlo Study of Quantized SU(2) Gauge Theory. Phys. Rev. 
D21, 2308 (1980). 

[90] K. Fabricius and O. Haan. Heat Bath Method for the Twisted Eguchi-Kawai 
Model. Phys. Lett. B143, 459 (1984). 

[91] N. Cabibbo and E. Marinari. A new Method for Updating SU(N) Matrices in 
Computer Simulations of Gauge Theories. Phys. Lett. B119, 387 (1982). 



Bibliography 



85 



[92] G. Bhanot. The Metropolis Algorithm. Rept. Prog. Phys. 51, 429 (1988). 
[93] J. Nishimura, 02 2002. Private Communication. 



